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SUMMARY 


In  theoretical  and  numerical  analyses  of  high-altitude  nuclear  plasmas 
and  ionospheric  nonnuclear  plasmas,  the  assumptions  are  often  made  that  the 
parallel  plasma  conductivity  is  infinite  and  that  the  electron  dynamics 
perpendicular  to  the  magnetic  field  are  governed  by  E  x  B  motion.  These 
assumptions  imply  that  perpendicular  electric  fields  map  perfectly  along  the 
magnetic  field  and  that  electrons  always  remain  in  the  same  flux  tube. 

Actually,  electron-ion  and  electron-neutral  collisions  assure  that  the 
classical  parallel  plasma  conductivity  is  finite,  albeit  generally  much  larger 
than  the  perpendicular  plasma  conductivity.  This  means  that  perpendicular 
electric  fields  do  not  map  perfectly  along  the  magnetic  field  and  that  electrons, 
which  are  originally  in  the  same  flux  tube,  do  not  always  have  to  remain  in  the 
same  flux  tube. 

The  primary  purpose  of  this  report  is  to  explicitly  demonstrate  that  the 
reality  of  finite  parallel  conductivity  can  substantially  impact  the  gross 
dynamics  of  high-altitude  plasma  clouds  as  well  as  the  structuring  of  the  clouds 
via  Rayleigh-Taylor  or  gradient-drift  mechanisms.  Explicit  quantitative 
demonstrations  are  given  for  the  three-dimensional  localization  of  electric 
field  patterns  around  plasma  clouds  and  for  the  fact  that  perpendicular  electric 
fields  do  not  necessarily  map  between  conjugate  points  along  the  geomagnetic 
field.  The  results  clearly  point  to  the  need  for  reconsidering  the  assumption  of 
infinite  parallel  conductivity  in  theoretical  and  numerical  models  of  the  high- 
altitude  nuclear  environment. 

A  secondary  purpose  of  this  report  is  to  illustrate  ways  that  ion-neutral 
diffusive  currents  and  ion-viscosity  can  affect  the  evolution  of  high-altitude 
plasmas. 

This  report  is  a  compilation  of  three  independent  chapters,  each  of  which 
demonstrates  that  finite  parallel  conductivity  and  associated  parallel  electron 
currents  can  directly  contribute  to  the  response  of  ionospheric  and  magneto- 
spheric  plasmas  to  an  applied  force.  The  third  chapter  should  be  of  special 
interest  to  the  reader  who  is  concerned  with  the  modeling  of  the  high-altitude 
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nuclear  environment,  as  the  chapter  includes  algebraic  algorithms  for  describing 
striations  in  the  high-altitude  nuclear  environment. 

Chapter  1,  "On  the  Response  of  Three-Dimensional  Plasma  Clouds  to  an 
Applied  Force,"  shows  how  a  three-dimensional  plasma  cloud  responds  to  an 
applied  current-density  field  resulting  from  an  applied  force.  Using  a  waterbag 
model  for  the  plasma  cloud,  it  is  demonstrated  that  parallel  electron  currents 
reduce  both  the  size  of  induced  perpendicular  electric  fields  and  associated 
E  x  B  velocities  when  the  applied  current-density  field  is  perpendicular  to  the 
ambient  magnetic  field.  It  is  suggested  that  anomalous  resistive  processes  can 
contribute  to  the  tendency  of  the  backside  of  plasma  clouds  to  steepen  even  if 
the  plasma  cloud  and  the  surrounding  plasma  are  Uniform.  Ion  viscosity  is 
explicitly  demonstrated  to  create  boundary  layers  near  the  edges  of  waterbag 
plasma  clouds.  These  boundary  layers  are  relatively  thin  at  low  altitudes  but 
btcome  increasingly  broad  at  higher  altitudes  if  ion-polarization  currents  are 
weak.  It  is  demonstrated  that  the  development  of  striations  can  be  hindered  when 
ion-neutral  diffusive  currents  contribute  to  the  overall,  nondivergent,  current- 
density  field.  Striation  evolution  is  discussed  in  terms  of  a  "U-shaped"  curve 
(actually  a  straight  line)  which  is  indicative  of  the  parameter  regime  when  ion- 
neutral  diffusive  currents  may  play  an  appreciable  role  in  the  dynamics  of 
striations  and  striation  "freezing."  A  parameter,  W,  is  defined  as  a  measure  of 
the  size  of  induced  electric  fields.  Ion-Pedersen,  ion-polarization,  and 
parallel-electron  resistive  currents  all  contribute  to  W. 

Chapter  2,  "The  Role  of  Finite  Parallel  Length  on  the  Onset  of 
Striations,"  uses  a  simple  plasma  model  to  show  that  the  finite  parallel  length 
of  ionospheric  plasma  clouds  plays  a  role  in  the  delayed  formation  of  striations 
by  tending  to  reduce  the  growth  rate  of  striation  instabilities.  Specifically, 
if  the  electron  density  in  the  plasma  cloud  is  much  larger  than  the  electron 
density  of  the  background  electrons,  the  length  of  the  plasma  cloud  enters 
through  the  field-line  integrated  electron  density  of  the  cloud.  Taking  into 
account  the  field-line  integrated  electron  density,  the  onset  of  appreciable 
striations  is  assumed  to  occur  when  the  growth  rate  of  the  striation 
instabilities  is  much  faster  than  the  rate  of  change  of  the  bulk  plasma.  This 
point  is  quantified.  The  finite  parallel  length  of  plasma  clouds  tends  to  favor 
the  growth  of  striations  with  finite  rather  than  with  zero  perpendicular  wave 
numbers. 


Chapter  3,  "Late-Time  Striation  Eigenmodes  Along  the  Geomagnetic  Field," 
derives  a  second-order  linear  differential  equation  for  determining  the  linear 
eigenmodes  of  striation  (i.e.,  gradient-drift  and  Rayleigh-Taylor)  instabilities 
along  the  geomagnetic  field.  Eigenvalues  (i.e.,  temporal  frequencies  in  complex 
space)  are  determined  along  with  the  eigenmodes.  The  differential  equation 
includes  the  contribution  of  two  significant  physical  effects:  ion-polarization 
currents,  and  inductive  electric  fields,  and  is  solved  for  parameters  and 
geometry  appropriate  to  late-time  HANE  plasmas,  as  described  by  MELT.  Ion- 
Pedersen,  ion-gravity,  and  ion-curvature  drifts  are  all  considered  as  possible 
contributors  to  the  formation  of  striations  at  late  times.  For  a  specified  late¬ 
time  plasma  configuration,  striations  tend  to  have  a  larger  scale  size,  perpen¬ 
dicular  to  the  geomagnetic  field,  at  higher  altitudes  than  at  lower  altitudes. 
The  ion-Pedersen  drift  becomes  an  especially  important  contributor  to  the 
formation  of  late-time  striations  with  small  perpendicular  scale  size. 
Electromagnetic  effects  are  shown  to  be  important  for  modes  with  smaller 
perpendicular  wave  numbers  (i.e.,  <1  km"^}  and  ion-viscosity  is  important  for 
modes  with  larger  perpendicular  wave  numbers  (i.e.,  ^  4  km-^).  Analytic 
criteria  are  given  for  the  perpendicular  wave  numbers  of  the  instability  at  which 
the  real  part  of  the  frequency  is  equal  to  the  growth  rate  and  at  which  there  is 
complete  stabilization.  It  is  argued  that  the  equating  of  the  real  part  of  the 
frequency  and  the  growth  rate  gives  a  meaningful  estimate  of  the  perpendicular 
wave  number  appropriate  to  striation  "freezing."  An  algorithm,  suitable  for 
implementation  in  modeling  codes,  is  given  for  determining  the  perpendicular 
wave  number  at  "freezing."  Simple  empirical  estimates  for  the  parallel  spatial 
extent  of  the  striation  density  fluctuations  are  also  described. 


PREFACE 
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Appendix  B  uses  (1.21a-b)  to  derive  the  total  electrostatic  potential 
for  the  case  when  the  external  applied  force  is  completely  electrostatic.  The 
appendix  demonstrates  that  the  magnitude  of  the  total  electric  field  inside  the 
ellipsoid  is  reduced  relative  to  the  value  outside  the  ellipsoid  if  the 
appropriate  conductivity  inside  the  ellipsoid  is  larger  than  the  corresponding 
conductivity  outside  the  ellipsoid.  For  example,  if  the  external  applied 
electric  field  is  purely  x-directed,  then  l°x</ox>|  >  1  is  required  in  order  that 
the  electric  field  inside  the  ellipsoid  be  smaller  in  magnitude  than  far  from  the 
ell ipsoid. 

The  depolarization  coefficients  [i.e.,  h  (0),  h  ( 0 ) ,  and  h  (0)]  satisfy 

A  Jr  £. 

[Landau  and  Lifshitz,  I960] 

0  <  hx(0)  <  hy(0)  <  hz(0)  ,  (1.24) 

if  a^  >  bt  >  ct  and  s  >  0.  Moreover  [Landau  and  Lifshitz,  I960], 

1  =  hx(0)  +  hy(0)  +  h2(0)  ,  (1.25) 


and  so 


0  <  hx(0,  hy(0,  hz(0  <  1  .  (1.26) 

It  is  clearly  evident  that  the  depolarization  factors  strongly  affect  the  size  of 
induced  electric  fields. 

In  general,  the  integrals  determining  the  depolarization  coefficients 
are  elliptic  and  of  the  second  kind  [Stratton,  1 94 1  ] .  However,  they  can  be 
evaluated  in  terms  of  more  elementary  functions  in  certain  limits.  Examples  are 
given  in  Table  1.2. 

The  first  two  entries  in  Table  1.2  correspond  to  the  ellipsoid  becoming  a 
long  elliptical  or  circular  cylinder  u-,  v-,  w-coordinate  system  defined  by 
(1.7).  The  entries  agree  with  previous  results  for  two-dimensional,  infinitely 
long  plasma  clouds  [Linson  and  Meltz,  1972].  These  entries  indicate  that  nonzero 
Joz£  do  not  9enerate  nonzero  electrostatic  potentials,  However,  nonzero  jQx> 


=  0,  (1.17)-(1.19)  give 


For  J'oy  =  J‘oz 


*<  =7 


x  (Jox<  ~  J'ox>)hx(0) 


u>  1  +(  — —  -  1)  ^  (o) 

°i> 


_  X  [_ J OX<  “  ^OX>^X^ 


°r'  1  +(=tl  -  1)  fix(0) 
°1> 


For  the  general  case  of  jQX,  j0  ,  and  jQZ  not  equal  to  zero,  symmetry 


indicate  that: 


x  (Jox<  •  Jox>)hx(0)  '■  (J—  -  )h-(0) 


(J)<  =  *  +  _JL  _ _  °y>  v 


Ji>  1  +[^i  -  1)  h  (0)  0i>  1  -  1)  h  (0) 

°1>  /  K>  '  " 


z  ^Joz<  "  ^oz>^hz^ 


z> 


1  + 


M 

az>  / 


h2(0) 


x  ^'ox<  ~  Jox>^x^  +  y  ^'oy<  ~  ’Vy^.y^ 


>  — 
a 


i> 


1  + 


M 


hx(0)  0l>  1  +(  ~  -  1J  hy(0) 

"1>  '  '  °l> 


+  _L_  (joz<  '  J'oz>^hz^^ 


z> 


(0) 


(1.20a) 


(1.20b) 


arguments 


(1.21a) 


(1.21b) 


(1.15) 


*0><  =  C^F(C)F(n)F(C) 


with 


c? 


-W°i>>1/2 


2  -2nl/2 


Kbt  ‘  at)(ct  ‘  at^ 


F(X)  =  (X  +  aj) 


1/2 


(1.16) 


The  electrostatic  potential  generated  by  the  ellipsoid  is 


=  C:  F(0F(n)F(C) 


C2  F(0F(n)F(C)hx(0 


(1.17a) 


(1.17b) 


with 


hx(0 


=  atbtct  f _ _ dr 

2  J  (r  +  a2)1,5(r  +  b? 


(r  +  aJ)A'3(r  +  b*)°'5(r  +  c2)0'5 


(1.18) 


At  the  surface  of  the  ellipsoid  (i.e.,  £  =  0),  the  electrostatic 
potential  and  the  normal  component  of  the  current  density  are  both  continuous 
[Stratton,  1941],  i.e.. 


4> 


5=0 


=  *< 


5=0 


(1.19a) 


K  + j_^o>\ 

\3^  ax>  /5- o 


^  !^  +  _l3_V 

^1>  3C  °i>  3C  /  5=0 


.  (1.19b) 


Equation  (1.19b)  uses  the  fact  that  the  appropriate  values  for  the  external, 
applied,  current-density  field  are  the  transformed  space  defined  by 
the  coordinates  u,  v,  and  w  [Landau  and  Lifshitz,  I960]. 


:: 
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2  2  2 
For  C  >  -  ct,  (1.12)  becomes  an  equation  for  an  ellipsoid.  For  -c£  >  n  >  -b£, 

2  2 

(1.12)  becomes  an  equation  for  a  hyperboloid  of  one  sheet.  For  -b£  >  G  >  -a£, 
(1.12)  becomes  an  equation  for  a  hyperboloid  of  two  sheets.  In  terms  of  £,  n, 
and  4,  the  cartesian  coordinates  u,  v,  and  w  are: 


(S  +  <Hn  *  aj;)(4  +  a‘) 
0>t  -  ‘t^t  -  at> 


(1.13a) 


v  =  ± 


(5  +  b|)(n  +  b^)(c  +  bj) 

7“2  Jw  2  TTT- 

(ct  -  bt)(at  -  bt) 


(1.13b) 


(5  *  c*)(n  +  c*)(S  +  c*) 
(a|  -  cf)(bf  -  cf)  _ 


(1.13c) 


It  is  assumed  that  the  external  applied  force  drives  a  uniform  current- 
density  field,  j  ,  both  far  from  the  ellipsoid  and  inside  the  ellipsoid 
[Stratton,  1941;  Landau  and  Lifshitz,  I960]: 


(1.14a) 


*0*  *  'W  -  -  W  •  (1'14b> 

In  (1.14a)  the  unit  vectors  in  the  x-,  y-,  and  z-directions  are  represented  by  x, 
y,  and  z,  respectively.  Laplace  transformation,  (1.5a),  permits  temporal 
variations  of  and  4^  to  be  taken  into  account. 

“►  /\ 

For  the  moment  the  special  case,  jQ^  =  jox^x  and  4>0^  =  -jox^x,  is 
considered.  This  simplification  is  subsequently  relaxed  using  symmetry 
arguments.  With  the  use  of  (1.5a)  and  (1.14a),  it  is  then  evident  that 
[Stratton,  1941] 
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In  the  transformed  coordinates  the  equation  determining  the  potential  for  the 
external  plasma  is  Laplace's  equation  which  is  appropriate  for  determining  the 
electrostatic  potential  in  a  vacuum.  Equation  (1.8b)  indicates  that  the 
conductivity  tensor  remains  a  diagonal  tensor  with  the  transformation  (1.7),  but 
the  elements  of  the  tensor  have  changed  (i.e.,  o  -*  a  /o¥.»  ou  ■*  a  /o  ,  and 

•k’VV1' 

The  surface  of  the  ellipsoid  is  determined  by  the  equation 
[Stratton,  1941]: 


(1.9) 


The  semiprincipal  axes  have  length  aQ,  bQ,  and  c^.  Using  (1.7)  permits  (1.9)  to 
be  rewritten  as 


Hence,  in  the  transformed  coordinates  the  plasma  cloud  is  still  an  ellipsoid  but 
the  semiprincipal  axes  are  now  at,  b^,  and  c^.. 

The  exact  functional  forms  for  ^  and  can  be  most  directly  determined 
using  the  ellipsoidal  coordinates  E,  n,  and  4  [Stratton,  1941;  Landau  and 
Lifshitz,  I960].  If  r  =  5,  n,  or  c ,  then  the  ellipsoidal  coordinates  each 
satisfy  the  equation 


(1.12) 
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ds  g  exp(st) 


(1.5b) 


9  =  J 

scr’“ 

In  (1.5a-b),  g  is  a  function  of  time  and  g  is  the  Laplace  transform.  The  contour 
indicated  by  (1.5b)  is  to  the  right  of  all  poles  of  g  in  the  complex  plane. 

If  <t>£  *  0  at  t  *  0,  then  the  application  of  (1.5a)  to  (1.4a-b)  results  in: 


0  = 


(1.6a) 


0  = 


+  a 

y< 


+ 


(1.6b) 


Appendix  A  describes  the  mathematical  forms  for  aY,  a„>  and  o,  used  in 

this  paper.  The  conductivities  in  the  Laplace-transform  space  (i.e.,  o¥,  ou,  and 
—  .  *  y 

az)  are  also  described  in  the  appendix. 

Consider  the  coordinate  transformation  on  page  62  of  Landau  and  Lifshitz 

[i960]  : 


u  =  x/(ax>)1/2  »  v  =  ^/(V)1/2  ’  W  =  z/(oz>)1/2  '  U-7)' 


With  (1.7),  (1.6a)  and  (1.6b)  can  be  rewritten  as 


The  electros  atic  fields  can  be  determined  from  the  potential,  $, 
through  the  relation, 

E  *  -  V  4>  •  (1.1) 

The  potential  itself  is  determined  from  the  constraint  that  the  divergence  of  the 
current  is  equal  to  zero,  i.e., 

0  =  7*o  •  V<|>  »  (1.2) 

where  o  is  the  conductivity  tensor.  Throughout  this  paper,  the  conductivity 
tensor  is  assumed  to  be  diagonal,  i.e., 


(1.3) 


Now  using  the  subscript  > (<)  to  denote  parameters  outside  (inside)  the 
ellipsoid  permits  (1.2)  to  be  written  as: 


(1.4a) 


0 


y< 


(1.4b) 


The  conductivities  and  potentials  in  (1.4a-b)  may  depend  on  time.  The 
temporal  variations  can  be  accounted  for  by  using  the  technique  of  Laplace 
transformation  [Spiegel,  1965]: 


g  =  J  dt  g  exp(-st) 
0 


(1.5a) 
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Figure  1.1.  Illustration  of  a  planar  section  of  an  ellipsoidal  plasma  cloud.  The 
section  is  parallel  to  the  x-y  plane.  The  magnetic  field  is 


-directed.  For  an  applied  y-directed  force. 


the  associated 


plasma  current  densities  are  The  induced  electric  field 
inside  the  cloud  points  in  the  negative  x-direction.  The  most 
likely  locations  for  anomalous  resistivity  and  striation  formation 
are  also  shown. 


Table  1.1.  Examples  of  external  perpendicular  forces  appropriate  to  uniform 
magnetic  fields.  Symbols  used  are:  mi  (ion  mass),  (perpendicular 
component  of  the  acceleration  of  gravity),  EQi  (perpendicular 
electric  field),  d/dt  (derivative  with  respect  to  time),  vin  (ion- 
neutral  collision  frequency),  and  v_  (neutral  velocity). 


"freezing."  From  the  mathematical  standpoint,  the  curve  is  obtained  by  equating 
the  real  and  imaginary  part  of  the  linear  frequency  associated  with  the  striation 
instabilities  [Sperling  and  Glassman,  1983;  Sperling,  1983a].  It  is  shown  that 
the  curve  obtained  in  this  way  is  consistent  with  the  observed  perpendicular 
scale  sizes  of  striations.  A  parameter,  W,  is  defined  as  a  measure  of  the  size  of 
induced  electric  fields.  Striation  evolution  is  discussed  from  the  viewpoint 
that  their  existence  must  be  indicative  of  the  elongation  and  motion  of  a  plasma- 
cloud  material  in  the  direction  of  an  applied  perpendicular  force. 

This  chapter  is  divided  in  the  following  way.  Using  Laplace  transforms, 
section  1.2  derives  the  induced  three-dimensional  electrostatic  potential 
resulting  from  the  application  of  an  applied  force  to  a  plasma  cloud  and  the 
surrounding  plasma.  Ion-Pedersen,  ion-polarization,  and  parallel  electron- 
resistive  currents  are  all  considered.  Section  1.3  discusses  the  inverse 
Laplace  transformation  of  certain  equations  in  .section  1.2.  Section  1.4 
discusses  applications  to  the  evolution  of  plasma  clouds  and  the  "U-shaped” 
curve.  Section  1.5  is  a  summary.  Six  appendices  clarify  the  discussion  in  the 
various  sections. 

1.2  CALCULATION  OF  INDUCED  ELECTROSTATIC  FIELDS 

A  relatively  simple  model  is  used  to  ascertain  the  role  of  three- 
dimensional  effects  on  the  electric  fields  generated  by  a  plasma  cloud  in  a 
uniform  current  field.  The  electric  fields  are  assumed  to  be  electrostatic  and 
the  current-density  field  is  assumed  to  be  the  result  of  an  external  force  (e.g., 
gravity,  drag  from  neutral  wind,  ion  inertia,  curvature  of  magnetic  field  lines, 
etc.)  imposed  on  the  plasma  cloud  and  the  surrounding  plasma.  Examples  of 
expressions  for  the  current-densities  due  to  external  forces  are  found  in 
Table  1.1,  page  50  of  Longmire  [1963],  and  page  673  of  Kilb  [1977]. 

The  plasma  cloud  is  assumed  to  have  the  general  configuration  of  an 
ellipsoid.  All  plane  sections  of  an  ellipsoid  are  ellipses  or  circles. 
Figure  1.1  is  an  illustration  of  a  planar  section  of  the  ellipsoid  parallel  to 
the  x-y  plane.  The  plasmas  in  the  interior  and  the  exterior  of  the  ellipsoidal 
cloud  are  assumed  to  both  be  uniform.  The  ambient  magnetic  field  is  uniform 
everywhere  and  z-directed. 
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Ion  viscosity  is  quantitatively  demonstrated  to  be  a  process  which  can 
affect  the  development  of  induced  electric  fields  within  plasma  clouds  when  the 
applied  force  acts  transverse  to  the  magnetic  field.  Viscosity  is  an  important 
consideration  because  the  E  x  B  velocitiec  associated  with  the  induced  electric 
fields,  can  be  nonuniform,  particularly  near  the  boundary  of  a  plasma  cloud. 
Specif ical ly,  it  is  demonstrated  that  viscous  boundary  layers  develop  near  the 
edge  of  a  waterbag  model  for  a  plasma  cloud.  Because  the  E  x  B  velocity  is 
nonuniform  within  the  boundary  layer,  the  boundary  layer  is  effectively  an 
ablation  layer  [Linson  and  Meltz,  1972].  For  steady-state  situations,  the 
boundary  layer  is  thinnest  at  low  altitudes  where  the  ion-neutral  collision 
frequency  is  large,  but  at  high  altitudes,  where  the  ion-neutral  collision 
frequency  is  small,  the  boundary  layer  is  broad.  Hence,  a  waterbag  model  may  not 
be  appropriate  for  high-altitude  plasma  clouds  near  steady  state.  However,  if 
ion-polarization  currents  are  significant,  thin  boundary  layers  can  form  at  high 
altitudes,  even  if  the  ion-neutral  collision  frequency  is  small. 

Ion-neutral  diffusion  can  reduce  the  tendency  for  induced  perpendicular 
electric  fields  to  form  in  the  ionosphere,  if  the  perpendicular  density  gradient 
is  such  that  the  ion-neutral  diffusive  current  plays  an  essential  role  in  the 
maintenance  of  overall,  nondivergent,  current-density  fields.  For  the 
appropriate  value  of  the  perpendicular  density  gradient,  it  is  not  necessary  to 
invoke  induced  electric  fields  and  consequently  the  presence  of  an  applied  force 
correspondingly  does  not  necessarily  result  in  overall,  E  x  B,  plasma  motion  and 
elongation  in  the  direction  of  the  applied  force.  Since  the  formation  of 
striations  is  a  consequence  of  the  tendency  of  a  plasma  cloud  to  move  and 
elongate  in  the  direction  of  a  force,  it  can  be  expected  that  the  development  of 
striations  is  hindered  by  ion-neutral  diffusion. 

Striation  evolution  has  been  discussed  in  terms  of  the  "U-shaped"  curve 
[Linson  ’nd  Meltz,  1972;  McDonald,  et  al.,  198l] .  Previous  analyses  have 
considered  the  curve  from  the  viewpoint  of  direct  balance  between  striation 
growth  and  ambipolar  diffusion,  but  this  balance  is  not  consistent  with  the 
observed  perpendicular  scale  sizes  for  striations  [McDonald,  et  al.,  198l] .  In 
this  chapter,  the  appropriate  curve  (actually  a  straight  line)  is  considered 
to  indicate  the  parameter  regime  when  ion-neutral  diffusive  currents  play  an 
appreciable  role  in  the  dynamics  of  striations  and  contribute  to  striation 
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1.  ON  THE  RESPONSE  OF  THREE-DIMENSIONAL  PLASMA 
CLOUDS  TO  AN  APPLIED  FORCE 


1.1  INTRODUCTION 

In  a  magnetized  plasma,  the  application  of  a  force  (e.g.,  gravity,  drag 
from  neutral  wind,  ion  inertia,  etc.)  transverse  to  the  ambient  magnetic  field 
results  in  a  current  density  transverse  to  both  the  ambient  magnetic  field  and 
the  force  [Longmire,  1963].  Considered  by  itself,  this  current  density  may 
diverge  as  it  flows  through  an  inhomogeneous  plasma.  Concomitant  charge 
accumulations  result  in  induced  electric  fields  transverse  to  the  ambient 
magnetic  field.  Plasma  motion  and  elongation  in  the  direction  of  the  force  are  a 
consequence  of  the  induced  electric  fields.  Indeed,  the  formation  of  striations 
is  one  way  that  plasma  clouds  tend  to  move  and  elongate  in  the  direction  of  a 
force  [Linson  and  Workman,  1970;  Volk  and  Haerendel,  1971;  Perkins,  et  al., 
1973]. 

This  chapter  examines  the  significant  contribution  of  parallel  electron 
currents  on  the  maintenance  of  overall  quasineutrality  when  an  applied  force  acts 
on  both  a  plasma  cloud  and  the  surrounding  plasma.  It  is  explicitly  demonstrated 
that  parallel  electron  currents  tend  to  short-circuit  the  electrostatic 
potential  when  the  applied  force  is  perpendicular  to  the  magnetic  field.  Hence, 
the  induced,  perpendicular  electric  fields,  within  a  plasma  cloud,  are  corres¬ 
pondingly  reduced  when  the  self-consistent  parallel  electron  currents  are  large. 

When  the  self-consistent  parallel  currents  are  so  large  that  the 
associated  electron  speed  is  larger  than  the  ion  thermal  or  electron  thermal 
speed,  turbulence  processes  or  double  layers  may  develop  and  tend  to  retard  the 
currents  [Block,  1972;  Papadopoulos ,  1977;  Alfven,  198l] .  Under  such 
circumstances,  the  parallel  electron  currents  are  less  effective  in  short- 
circuiting  induced  electrostatic  potentials.  It  is  demonstrated  that  the 
presence  of  anomalous  resistive  processes  can  contribute  to  the  tendency  of  the 
backside  of  ionospheric  barium  clouds  to  steepen  even  if  the  plasma  cloud  and 
the  surrounding  plasma  are  initially  uniform.  The  tendency  of  anomalous 
resistive  processes  to  be  important  increases  with  altitude. 


and  j'0y^  may  generate  nonzero  electrostatic  potentials,  <J>^,  depending  on  the 
precise  geometry  of  the  elliptical  cylinder.  For  example,  consider  t>t  »  at  then 
the  elliptic  cylinder  becomes  a  thin  sheet  with  normal  in  the  x-direction.  In 
this  case,  j.v>  generates  a  substantial  electrostatic  field  since  h  (0)  =  1. 
However,  hy(0)  =  0  and  so  j  ^  does  not  generate  a  substantial  electrostatic 
field.  Evidently,  a  relatively  large  surface  area  normal  to  the  applied  current 
is  required  for  generating  an  appreciable  electrostatic  field. 

In  general,  it  can  be  expected  that  the  finite  extent  of  plasma  clouds 
along  the  magnetic  field  should  have  a  meaningful  effect  on  the  generation  of 
perpendicular  electrostatic  potentials  if  ct  <  at,b^.  Under  such  circum¬ 
stances,  the  ellipsoid  in  the  u-,  v-,  and  w-coordinate  system  has  a  relatively 
short  extent  along  the  magnetic  field.  For  ct  «  at,bt. 


dr 


(r  ♦  c2)3'2 


=  1 


(1.27) 


The  combination  of  (1.25)  and  (1.27)  gives 

hx(0)  =  hy(0)  =  0  (1.28) 

Equations  (1.27)  and  (1.28)  are  in  agreement  with  entry  3  of  Table  1.2  if 
ct  <<:  at,bt’  **  then  f°llows  if  ct  <<:  at’bt’  an  aPPreciable  electrostatic 

potential  is  not  generated  by  nonzero  Jox^>  J0y^'  This  t)ecause  parallel 
electron  currents  play  a  central  role  in  maintaining  overall  quasineutrality  and 
tend  to  short-circuit  the  potential  differences  which  would  arise  if  perpen¬ 
dicular  currents  alone  were  important.  However,  the  presence  of  nonzero  jQZ^  may 
generate  an  electrostatic  field. 

In  the  limit  when  two  of  the  semiprincipal  axes  are  equal  (i.e.,  at  =  bt, 
at  =  ct’  or  bt  =  ct^  tbe  eliiPSoids  become  spheroids.  Specifically,  if  c^  <  at  = 
bt  the  plasma  cloud  is  an  oblate  spheroid  but  if  ct  >  a^  *  b^  the  plasma  cloud  is 
a  prolate  spheroid. 

Figure  1.2  is  a  plot  of  ^x(0)  and  Fz(0)  for  the  case  of  at  =  bt  and 
entries  4  and  5  in  Table  1.2  [Stratton,  1941;  Landau  and  Lifshitz,  I960].  In 


23 


agreement  with  entry  2  of  Table  1.2,  hv(0)  ■*  0.5  and  h  (0)  ■*  0  as  c*/a++®. 

p  p  _  _  Z  t  l 

Similarly,  if  c^/a^  -*■  0  then  hx ( 0)  -*■  0  and  hz(0)  •+  1  in  agreement  with  (1.27), 
(1.28),  and  entry  3  of  Table  1.2.  For  ct  =  at  -  b^  the  spheroid  becomes  a  sphere 
and  entry  6  in  Table  1.2  follows  directly  from  symmetry  arguments. 

2  7  2 

Far  from  the  ellipsoid  (i.e.,  £  »  a£,  b£,  c£). 


2  2  2  2 

P6  =  u^  +  v^  +  w^  =  - — +  1 


(1.29) 


hx(«  «  Fy(d  -  h2(5)  -  -ffc 


atbtCt 


(1.30) 


The  combination  of  (1.21b),  (1.29),  and  (1.30)  gives 


f1 


hx(0) 


1  +  -r-  -  1  MO) 


^Jox<  ~  Jox>^ 


(1.31) 


It  is  evident  from  (1.31)  that  the  electrostatic  potential  has  a  dipolar 
character  sufficiently  far  from  the  ellipsoidal  cloud  [Landau  and  Lifshitz, 
I960],  The  P"^  dependence  indicates  a  corresponding  rapid  spatial  decay  of  the 
electric  field. 


1.3  INVERSE  LAPLACE  TRANSFORMATION 

To  determine  the  electrostatic  potential  in  real  space,  the  inverse 
Laplace  transformation,  (1.5b),  can  be  applied  to  (1.21a-b).  In  general,  an 
analytic  expression  for  the  inverse  transformation  cannot  be  obtained  except  in 
certain  special  limits.  The  special  limits  emphasized  in  this  section  are 

ct  >>  at’bt  and  ct  <<c  at  =  br  11  is  assumed  th»t  Joy^=  =  0.  It  is  also 
assumed  that  is  constant  in  time  so  that 


Jox^  ~  Jox£  s 


(1.32) 


With  (1.21a),  (1.23a),  (1.32),  and  (A4)  of  Appendix  A,  the  final- 
value  theorem  for  Laplace  transformations  gives  [Spiegel,  1965] 


<t>  (t  <*)  =  lim  (s<J>  )  =  - — 

_  n  <'  Qt  V.  +  (ot  V- 
s-*0  >  in>  '  <  in< 


X(jOX<  “  ^OX>]bX^^ 


-  a>vin>)hx(°) 


,  (1.33a) 


j  (t  -*■  «o)  s  i  jm  -  (jox>  “  jox<)ot<vin<hxt0)  n 

X<(  *  .  "o  (SJ<)  Vln>  +  lVin<  -  Vin>)M°>  ’ 


hx(0) 


a  v. 

>  in> 


a  v. 

>  in> 


(1.34) 
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If  ct  »  at,bt,  then  the  ellipsoid  becomes  a  very  long  elliptical  cylinder  in  the 
transformed  coordinates  specified  by  (1.7).  Parallel  electron  currents  do  not 
play  a  major  role  in  maintaining  overall  quasineutrality.  The  combination  of 
(1.21a),  (1.32),  (A4)  of  Appendix  A,  and  entry  1  of  Table  1.2  gives 
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with  determined  from  (A5)  of  Appendix  A.  The  combination  of  (1.5b)  and 
(1.34)  gives 
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It  follows  from  (1.36)  that 
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If  ct  «  a^.  =  bt,  then  the  ellipsoid  becomes  a  thin  circular  disk  in  the 
transformed  coordinates  specified  by  (1.7).  Parallel  electron  currents  play  an 
essential  role  in  maintaining  overall  quasineutrality  and  impede  the  tendency 
for  substantial  electrostatic  potentials  to  form.  The  combination  of  (1.21a), 
(1.31),  (A4)  of  Appendix  A,  and  entry  4  of  Table  1.2  gives 
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Equation  (1.38)  may  be  rewritten  as 
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The  combination  of  (1.5b),  (1.32),  and  (1.40)  gives  [Spiegel,  1965] 
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Furthermore,  if  |x+j  «  vin>  and  a</a>  >  1  then  <  0.  It  then  follows  that  <\>< 
has  the  following  form  for  t  »  v.n>: 
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When  neutrals  are  sparse  and  becomes  very  small,  then  1^1  can  tend  to 
become  large.  However,  a  large  value  for  the  parallel  conductivity  counteracts 
the  formation  of  large  electrostatic  potentials  by  providing  a  mechanism  for 
short-circuiting. 

For  ct  «  at  =  b^,  the  maximum  magnitude  for  the  steady-state  parallel 
electron  current  on  the  surface  of  the  ellipsoid  is  determined  from  Appendix  C 
to  be 
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This  expression  assumes  that  =  joz^  *  0.  The  maximum  parallel  electron 

speed  corresponding  to  (1.45)  is 
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It  is  evident  from  (1.46)  that  |v  |  diverges  for  a,.  -*•«>  and  v.  -*■  0. 

If  ct  «  at,bt  then  parallel  electron  currents  are  important  and  tend  to 
reduce  the  effective  value  of  .  However,  it  follows  that  if  parallel 
electron  currents  are  an  essential  component  of  the  total  current  system, 
including  the  plasma  cloud,  then  the  parallel  currents  may  be  sufficiently  large 
to  generate  turbulence  or  double  layers  [Block,  1972;  Papadopoulos,  1977; 
Alfven,  1981].  Since  parallel  conductivity  is  effectively  decreased  by 
turbulence  and  double  layers,  it  follows  that  parallel  currents  are  corres¬ 
pondingly  decreased  and  less  able  to  contribute  to  the  overall  current  systems. 
An  increase  in  the  magnitude  of  |<j>>|  is  also  implied. 


If  the  effective  parallel  conductivity  is  reduced  by  plasma  processes 
(e.gl,  turbulence  or  double  layers),  (1.46)  correctly  demonstrates  that  there  is 
a  corresponding  reduction  in  the  parallel  electron  velocity,  the  free-energy 
source  responsible  for  the  turbulence  and  double  layers.  Equations  (1.45)  and 
(1.46)  also  show  that  small  values  for  the  magnitude  of  the  perpendicular 
conductivity  tend  to  increase  the  likelihood  of  large  parallel  electron  current 
densities  and  velocities.  This  is  a  consequence  of  perpendicular  current 
densities  being  unable  to  maintain  nondivergent  current-density  fields  by 
themselves.  Hence,  parallel  currents  must  contribute  in  a  substantial  way  to  the 
maintenance  of  overall  quasineutrality.  The  tendency  for  large  parallel 
electron  velocities  to  occur  is  clearly  increased  if  |j  |  »  |j  |  and  if  the 
electron  density  is  much  larger  inside  the  plasma  cloud  than  outside  the  cloud. 

Up  to  now  it  has  been  assumed  that  perpendicular  current  densities, 
associated  with  a  plasma  cloud  and  divergent  jQ,  are  of  the  ion-Pedersen  or  ion- 
polarization  types.  However,  ion  viscosity  (cf.  Appendix  D)  and  ion-neutral 
diffusion  (cf.  Appendix  E)  may  also  contribute  to  the  maintenance  of  overall 
current-density  fields  which  are  nondivergent.  Appendix  D  points  out  that  ion 
viscosity  permits  large  perpendicular  gradients,  at  lower  altitudes  in  the 
ionosphere  where  ion-neutral  collisions  occur  at  a  relatively  rapid  rate.  At 
higher  altitudes  where  ion-neutral  collisions  are  negligible,  ion  viscosity 
still  permits  the  formation  of  sharp  perpendicular  density  gradients  if  ion- 
polarization  currents  are  sufficiently  strong.  Appendix  E  shows  that  a 
particular  value  for  the  perpendicular  density  profile  within  a  one-dimensional 
plasma  cloud  permits  ion-neutral  diffusive  currents  to  contribute  to  overall 
nondivergent  current-density  fields  without  requiring  the  formation  of  electric 
fields.  Appendix  E  also  shows  that  the  required  scale  size  of  the 
perpendicular  density  gradient  is  consistent  with  the  observed  scale  size  of 
"frozen"  striations. 

1.4  APPLICATION  TO  THE  EVOLUTION  OF  PLASMA  CLOUDS  AND  THE  "U-SHAPED"  CURVE 

With  the  application  of  a  force,  a  plasma  cloud  responds  in  two  ways. 
First,  the  cloud  tends  to  polarize  in  such  a  way  that  the  motion  of  the  cloud 
minimizes  the  force  in  the  cloud  frame  of  reference.  For  example,  in  the 
ionosphere,  a  plasma  cloud  would  try  to  move  at  the  neutral  velocity  as  a  result 


of  ion-neutral  drag.  Analogously,  in  a  nondissipative  medium,  a  plasma  cloud 
would  tend  to  accelerate  in  accordance  with  the  gravitational  force.  Second, 
plasma  clouds  deform  through  elongation  in  the  direction  of  an  applied  force.  In 
this  way  the  surface  area  projecting  in  the  direction  of  a  force  is  minimized  and 
the  cloud  becomes  aligned  in  the  direction  of  the  applied  force. 

Figure  1.1  illustrates  -a  planar  section  of  a  three-dimensional 
ellipsoidal  plasma  cloud  parallel  to  the  z  =  0  plane.  Assume  that  the  electron 
density  within  the  plasma  cloud  is  constant  and  that  a  force  is  applied  in  the 
y-direction.  As  shown  in  Figure  1.1,  the  associated  current-density  field  is  in 
the  x-direction  for  ambient  magnetic  field  pointing  out  of  the  paper.  The 
induced  electric  field  within  the  plasma  cloud  points  in  the  negative  x-direction 
and  so  the  reasonable  result  is  obtained  that  the  associated  E  x  B  velocity  and 
the  applied  force  both  point  in  the  positive  y-direction. 

If  the  volume  of  the  plasma  cloud  remains  fixed,  it  is  evident  from  the 
discussion  in  section  1.3  and  the  entries  in  Table  1.2,  that  an  increase  in  the 
ratio,  bt/at,  implies  elongation  in  the  direction  of  an  applied,  y-directed 
force.  Moreover,  larger  values  for  the  ratios,  ct/at  and  ct/bt,  indicate  a 
reduction  in  the  short-circuiting  involving  parallel  electron  currents  and  also 
indicate  that  a  larger  surface  area  projects  in  the  direction  of  the  current- 
density  field  associated  with  an  applied,  y-directed  force.  Elongation  of  a 
plasma  cloud  along  the  magnetic  field  or  anomalous  resistivity  both  can 
contribute  to  larger  values  for  c^.  As  b^.  and  c^  increase  relative  to  a^,  larger 
values  for  hx(0)  are  implied  and  so  both  the  induced  electric  field  inside  the 
plasma  cloud  as  well  as  the  associated  y-directed  E  x  B  velocity  both  corres¬ 
pondingly  increase.  If  the  volume  of  the  plasma  cloud  remains  fixed  and  the 
electron  density  remains  uniform,  continued  acceleration  of  the  plasma  cloud  in 
the  direction  of  an  applied  force  can  be  expected  to  continue  until  hx(0)  =  1 
(i.e.,  ct  »  bt  »  at).  When  hx(0)  =s  1  further  elongation  and  acceleration  of 
material  in  a  plasma  cloud  can  occur  only  if  local  density  enhancements  develop 
which  effectively  increase  the  local  value  of  the  perpendicular  conductivity 
within  the  cloud  material.  Induced  electric  fields  should  be  correspondingly 
larger  at  the  locations  of  these  density  enhancements. 

In  order  for  the  overall  current-density  field  to  be  nondivergent,  it  is 
immediately  evident  from  Figure  1.1  that  the  induced  parallel  electron  current 


density  must  be  largest  near  the  two  points  where  the  surface  of  the  ellipsoid 
intersects  the  x-axis.  For  c£/a£  «  1  [e.g.,  Cg  =  10  km,  ag  =  2  km,  and 
(o2>/a>v^n>)  =  1  x  106],  (1-45)  provides  an  estimate  for  the  magnitude  of  this 
parallel  electron  current  and  (1.46)  provides  the  corresponding  estimate  for  the 
parallel  electron  speed.  For  representative  ionospheric  conditions  at  200  km 
and  400  km  altitude  as  illustrated  by  Goldman  et  al.  [1976],  Table  1.3 
quantifies  (1.46).  The  table  shows  that  the  parallel  electron  speed  can  exceed 
the  ion  thermal  speed  both  at  200  km  and  400  km  altitude.  In  fact,  the  parallel 
electron  speed  is  especially  large  at  an  altitude  of  200  km  where  it  may  even 
exceed  the  electron  thermal  velocity.  This  suggests  that  anomalous  resistive 
effects  (e.g.,  turbulence  and  double  layers)  [Block,  1972;  Papadopoulos,  1977; 
Alfven,  198l]  can  play  an  important  role  in  the  dynamics  of  the  plasma  cloud  by 
effectively  decreasing  the  parallel  electron  conductivity  of  the  external  plasma 
near  the  locations  where  the  plasma  cloud  intersects  the  x-axis  in  Figure  1.1. 
At  these  locations,  parallel  electron  current  densities  are  less  able  to 
contribute  to  overall  nondivergent  current-density  fields  than  locations  away 
from  the  x-axis.  (By  contrast  the  discussion  of  sections  1.2  and  1.3  assumed 
that  the  parallel  conductivity  is  the  same  everywhere  outside  the  plasma  cloud.) 
Consequently,  inside  the  cloud,  the  induced  electric  fields  and  the  associated 
induced  current  densities  must  be  larger  near  y  *  0  than  away  from  y  =0.  The 
attendant  E  x  B  velocity  is  also  largest  near  y  =  0.  It  follows  that  the  density 
of  the  cloud  cannot  remain  constant  but  must  tend  to  become  larger  for  y  >  0 
rather  than  y  <  0.  There  is  steepening  of  the  density  profile  for  y  >  0,  the 
backside  of  the  cloud,  and  the  density  profile  becomes  less  steep  for  y  <  0,  the 
frontside  of  the  cloud.  In  short,  turbulence  and  double  layers  can  contribute 
to  observed  steepening  and  accumulations  of  density  on  the  backside  of  plasma 
clouds  (and  also  striations)  by  permitting  inhomogeneous  E  x  B  velocities  within 
the  plasma  cloud.  Localized  anomalous  resistivity  is  a  mechanism  for  breaking 
symmetry  in  a  plasma  cloud  and  permitting  an  initially  uniform  plasma  cloud  to 
become  nonuniform. 

The  formation  of  striations  is  the  oft  observed  way  that  plasma  clouds 
tend  to  accelerate  and  elongate  in  the  direction  of  an  applied  force.  The 
approximate  location  where  striations  form  is  illustrated  in  Figure  1.1.  When 
striations  form,  it  can  be  expected  that  their  perpendicular  surfaces  would  tend 
to  be  tangent  to  an  applied  perpendicular  force  so  that  the  surface  area 
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Table  1.3.  |v_  as  determined  from  (1.46)  for  altitudes  of  200  and  400  km. 

Z  ^  j  HiQa 

It  is  assumed  that  the  applied  force  is  electrostatic  and  that  cE  /B 

4  ox 

=  1  x  10  cm/s.  All  perpendicular  plasma  currents  are  of  the  ion- 
Pedersen  type.  Plasma  and  neutral  parameters  are  determined  from 
Figures  2-4  of  Goldman  et  al.  [1976]  for  solar  minimum  and  solar 
maximum  conditions.  For  an  ion  atomic  mass  of  16  and  T.  =  T  = 

4  '  ^ 

0.1  eV,  the  ion  and  electron  thermal  speeds  are  7.7  x  10  cm/s  and 

7  3 

1.3  x  10  cm/s,  respectively.  It  is  assumed  that  n</n>  is  1  x  10 

and  10  for  altitudes  of  200  km  and  400  km,  respectively. 


Altitude  (km) 


lvz>,max^  (cm/s) 

1.6  x  108 

8.5  x  10 

solar  minimum 

'vz>,n,axl  <CI”/S> 

1.5  x  108 

1.6  x  10 

solar  maximum 

'  *1 


projecting  in  the  direction  of  the  force  is  minimized  (see  Figure  1.3). 
Striations  tend  to  be  aligned  in  the  direction  of  the  applied  force.  A 
completely  consistent  picture  is  obtained  if  the  induced  electric  fields  tend  to 
be  normal  to  the  applied  force  and  to  the  perpendicular  surfaces  of  the 
striations.  With  this  configuration  bulk  plasma  motion  occurs  in  the  direction 
of  the  applied  force  in  agreement  with  Prettie  et  al.  [  1983]  and  not  in  the 
direction  of  the  induced  electric  fields.  Hence,  striations  cannot  easily 
separate  from  one  another  and  they  appear  to  move  in  unison  in  the  E  x  B 
direction  [McDonald  et  al.,  1981 ] .  If  the  direction  of  the  applied  force  should 
change,  from  the  initial  form  indicated  in  Figure  1.3,  then  the  plasma  would  tend 
to  change  from  the  configuration  shown  in  Figure  1.3  to  a  new  configuration  which 
is  aligned  in  the  new  direction  of  the  applied  force.  New  striations  may  develop 
as  the  new  configuration  is  attained. 

The  purpose  of  the  rest  of  this  section  is  to  explicitly  demonstrate 
additional  ways  that  three-dimensional  effects  contribute  to  the  evolution  of 
plasma  clouds  and  the  interpretation  of  the  evolution  using  the  "U-shaped"  curve 
[Linson  and  Meltz,  1972;  McDonald  et  al,  1981].  Steady-state  and  temporally 
varying  situations  are  discussed  consecutively.  For  simplicity,  it  is  assumed 

that  VS  *  -W* 

a.  Steady  State 

If  j  is  a  result  of  ion-Pedersen  drifts,  then 
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The  combination  of  (1.33b)  and  (1.47)  gives  the  total  x-directed  plasma  current 
within  the  ellipsoidal  plasma  cloud: 
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with  Mq  being  the  ratio  of  local  plasma  conductivities,  i.e. 
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The  geometry  of  the  p.asma  cloud  contributes  to  (1.48)  through  the  depolarization 
factor,  hx(0). 

Two-dimensional  analyses  of  striation  evolution  in  barium  clouds  often 
use  field-line  integration  of  Pedersen  conductivity  to  take  the  distribution  of 
plasma,  along  the  geomagnetic  field,  into  account  [Linson  and  Meltz,  1972; 
McDonald  et  al.,  1981;  Prettie  et  a!.,  1983].  Appendix  F  shows  that  the 
concept  of  field-line  integration  is  appropriate,  to  within  a  factor  of  order 
unity,  when  ct,  a^  «  bt  or  c^  «  at  =  b^.  but  is  not  appropriate  to  c t  ~  at  or 


Now  two-dimensional  studies  and  entry  2  of  Table  1.2  give  [Linson  and 
Meltz,  1972;  Overman  et  al.,  1983] 


Jx<,tot 


(1.50) 


for  an  infinitely  long  circular  cylinder.  By  analogy  define  a  parameter  W  such 


Jxc,tot 
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It  then  follows  that 


W  =  1  +  2(Mq  -  l)hx(0) 


(1.52) 
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the  case  in  the  ionosphere  [Goldman  and  Sperling,  1982].  Third,  parallel 
electric  fields  can  be  mainly  inductive,  mainly  electrostatic,  or  both  inductive 
and  electrostatic.  The  inductive  contribution  to  the  electric  fields  implies 
modal  magnetic  fields  orthogonal  to  the  ambient  magnetic  field. 

Artificial  plasma  clouds  in  the  ionosphere  have  a  finite  spatial  extent 
that  can  influence  the  development  and  properties  of  plasma  instabilities.  In 
this  chapter,  particular  emphasis  is  placed  on  determining  certain  effects  on  the 
Rayleigh-Taylor  and  E  x  B-gradient-drift  class  of  instabilities  resulting  from 
the  finite  length  of  plasma  clouds  along  the  geomagnetic  field.  The  theory  is 
applied  to  the  onset  of  striations  in  barium  clouds. 

It  is  demonstrated,  with  a  rather  simple  plasma  geometry,  that  finite 
plasma  length  implies  parallel  modal  currents  and  parallel  modal  electric  fields 
which  contribute  to  the  development  of  eigenmode  structure  along  the  geomagnetic 
field.  Specifically,  the  finite  field-line  integrated  electron  density  of 
ionospheric  plasma  clouds  plays  a  role  in  the  delayed  formation  of  striations  by 
tending  to  reduce  the  growth  rate  of  striation  instabilities.  Taking  into 
account  the  field-line  integrated  electron  density,  the  onset  of  appreciable 
striations  is  assumed  to  occur  when  the  growth  rate  of  the  striation  insta¬ 
bilities  is  much  faster  than  the  rate  of  change  of  the  bulk  plasma.  This  point  is 
subsequently  quantified.  The  finite  parallel  length  of  ionospheric  plasma 
clouds  tends  to  favor  the  growth  of  striations  with  finite  rather  than  with  zero 
perpendicular  wave  numbers. 

In  section  2.2,  the  calculation  model  and  its  implications  are  detailed. 
Section  2.3  describes  the  quantitative  evaluation  of  the  model  for  parameters 
appropriate  to  barium  releases.  Concluding  remarks  are  given  in  section  2.4. 

2.2  CALCULATION  MODEL 

The  plasma  model  is  very  similar  to  the  one  used  in  Sperling  [l983c]  to 
study  the  effect  of  finite  parallel  length  on  the  loss-cone  instability 
hypothesized  to  be  responsible  for  prompt  striations.  The  finite  extent  of  the 
unstable  plasma  cloud  is  approximated  by  assuming  that  it  is  localized  to 
\z\  <  Zq  (see  Figure  2.1).  Here  z  is  the  coordinate  along  the  constant  magnetic 
field,  "b,  and  Zq  is  a  positive  constant.  The  plasma  for  |z|  <  Zq  has  a 
perpendicular  density-gradient  scale  length,  Lj_,  but  the  plasma  for  |z|  >  Zq  is 
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2.  THE  ROLE  OF  FINITE  PARALLEL  LENGTH  ON  THE  ONSET  OF  STRIATIONS 


2.1  INTRODUCTION 

In  the  disturbed  ionosphere,  Rayleigh-Taylor  [Scannapieco  and  Ossakow, 
1976;  Ott,  1978;  Ossakow  et  al.,  1979;  Zalesak  and  Ossakow,  1980]  and  E  x  B 
gradient  drift  instabilities  [Simon,  1963;  Linson  and  Workman,  1970;  McDonald 
et  al.,  1980,  198l]  are  considered  to  be  primary  sources  for  fluid  structuring. 
These  instabilities  are  predicted  to  be  aligned  along  magnetic  field  lines  with 
the  result  that  analysis  and  numerical  simulation  has  emphasized  the  flute 
approximation,  which  neglects  the  explicit  dependence  of  modes  on  the  coordinate 
parallel  to  the  ambient  magnetic  field.  Plasma  and  magnetic  field  variations  in 
the  direction  of  the  magnetic  field,  when  considered,  have  been  taken  into 
account,  in  a  gross  sense,  by  averaging  plasma  parameters  over  magnetic  field 
lines,  but  without  considering  the  implications  on  the  ionosphere  of 
accompanying  modal  variation  along  magnetic  field  lines.  One  exception  to  this 
generality  is  the  paper  by  Goldman  et  al.  [1976],  which  calculates  eigenmodes  in 
the  electrostatic  approximation,  ■  recognizing  that  modes  must  vary  along  the 
magnetic  field  lines  as  one  moves  away  from  the  source  of  instability,  an 
artificial  plasma  cloud. 

The  zero  electron  mass  approximation  is  often  made  in  the  theoretical 
analysis  of  plasma  oscillations  for  frequencies  below  the  lower-hybrid  frequency 
because  it  is  both  convenient  and  simplifying.  In  combination  with  the 
electrostatic  approximation,  the  zero  electron-mass  approximation  implies  that 
magnetic  field  lines  are  equipotentials  and  that  modes  have  a  flutelike 
character.  However,  the  approximation  of  zero  electron  mass  is  inconsistent  with 
the  observation  of  parallel  electric  fields  in  many  plasmas,  and  precludes 
consideration  of  a  number  of  associated  possible  effects.  First,  although  often 
much  smaller  than  other  electric  field  components,  parallel  electric  fields  can 
contribute  in  a  most  significant  way  to  the  properties  and  structure  of 
eigenmodes  [Allis  et  al.,  1963;  Paoloni,  1975;  Sperling,  1982a].  Second, 
parallel  electric  fields  can  result  in  the  transfer  of  energy  and  momentum 
between  the  mode  and  parallel  electron  or  ion  motion.  Depending  on  circum¬ 
stances,  this  coupling  can  be  collisionless  [Stix,  1962]  or  collisional,  as  often 
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turbulence  or  double  layers  may  develop  [Block,  1972;  Papadopoulos,  1977; 
Alfven,  198l] .  As  a  result,  the  tendency  of  parallel  electron  currents  to 
contribute  overall  current  patterns  can  then  be  reduced.  It  is  suggested  that 
this  anomalous  resistivity  contributes  to  the  breaking  of  symmetry,  even  in  the 
idealized  plasma  model  considered  in  section  1.2.  As  a  result,  the  anomalous 
resistivity  contributes  to  the  steepening  of  the  backside  of  plasma  clouds. 

Ion  viscosity  and  ion-neutral  diffusion  can  effect  the  formation  of 
induced  electric  fields  when  an  applied  perpendicular  force  acts  on  both  a  plasma 
cloud  and  the  surrounding  plasma.  For  steady-state  situations  and  when  the  ion- 
neutral  collision  frequency  is  large,  viscous  boundary  layers  develop  near  the 
edge  of  the  plasma  cloud.  For  a  barium  cloud  at  an  altitude  of  200  km,  the 

3 

boundary  layer  can  be  approximately  2  x  10  cm  thick.  Ion-neutral  diffusive 
currents  can  reduce  the  size  of  the  induced  electric  fields  required  to  maintain 
overall  nondivergent  current-density  fields.  Indeed,  no  induced  electric  fields 
are  required  if  the  density  profile  inside  a  plasma  cloud  has  the  form  indicated 
in  Appendix  E. 

Striation  evolution  has  been  discussed  in  terms  of  the  "U-shaped"  curve 
[Linson  and  Meltz,  1972;  McDonald  et  al.,  198l] .  In  the  present  chapter,  the 
curve  (actually  a  straight  line)  is  considered  to  indicate  the  parameter  regime 
when  ion-neutral  diffusive  currents  can  play  an  appreciable  role  in  the  dynamics 
of  striations  and  contribute  to  striation  "freezing"  [Sperling  and  Glassman, 
1983;  Sperling,  1983a].  A  parameter,  W,  is  defined  as  a  measure  of  the  size  of 
induced  electric  field.  Striation  evolution  is  discussed  from  the  premise  that 
their  existence  must  represent  the  means  whereby  a  plasma  cloud  elongates  and 
moves  in  the  direction  of  an  applied  force. 
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the  combination  of  ( 1 . 62 ) -( 1 . 64 )  gives 
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If  (1.47)  is  valid  then  -  MqA.  It  is  evident  from  (1.65)  that  steady-state 


values  for  W  are  attained  after  several  ion-neutral  collision  times  (i.e.,  v.  t 

in> 


and  vin<t  *  !)• 


As  altitude  increases,  the  ion-neutral  collision  frequency  decreases  and 
so  discrete  plasma  structures  should  be  able  to  accelerate  for  longer  periods  of 
time  at  higher  altitudes  than  at  lower  altitudes  where  ion-neutral  collisions  can 
rapidly  retard  plasma  motion  [Scholer,  1970;  Sperling,  1983b],  Indeed,  during  a 
magnetospheric  barium  release  striations  were  observed  to  accelerate  over 
periods  of  several  minutes  [Mende,  1973].  It  was  suggested  by  Mende  [1973]  that 
electromagnetic  effects,  not  considered  in  this  chapter,  can  contribute  to 
the  eventual  breaking  of  the  acceleration  at  high  altitude. 


1 . 5  SUMMARY 

A  very  general,  three-dimensional  ellipsoidal  model  for  plasma  clouds 
has  been  used  to  explicitly  determine  the  significant  role  of  parallel 
conductivity  and  self-consistent  parallel  electron  currents  on  the  formation  of 
induced  electric  fields  when  an  applied  force  acts  upon  the  plasma  cloud  and  the 
surrounding  plasma.  Field-line  integration  is  not  required  for  the  analysis  but 
the  results  (see  Appendix  F)  do  indicate  the  parameter  regime  where  field-line 
integration  is  approximately  appropriate. 


Self-consistent  parallel  electron  currents  in  the  three-dimensional 
configuration  may  tend  to  be  so  large  that  the  associated  parallel  electron 
velocities  may  exceed  the  ion-thermal  or  even  the  electron-thermal  speed  and  so 
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Both  ion-Pedersen  and  ion-polarization  currents  contribute  to  (1.59).  The  total 
x-directed  current  density  in  the  plasma  cloud  is 


^x<,tot  ~  ^ox<  +  ^x< 


If  in  analogy  to  (1.51)  a  parameter  W  is  defined  such  that 
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then  the  combination  of  (1.59)  and  (1.60)  gives 
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By  permitting  temporal  variations,  the  parameter  W  has  been  given  a  more  general 
form  than  indicated  by  (1.52). 

The  limit  of  two-dimensional  plasma  clouds  (i.e.,  ct  »  at,bt)  yields 
the  simplest  example  for  the  application  of  (1.62).  Equations  (1.36)  and  (1.59) 
then  give 
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Equation  (1.63)  implies  (1.48)  in  the  limit  of  t  -*•  ®  and  hx(0)  =  bg/(ag  +  bQ)  as 
suggested  by  entry  1  of  Table  1.2.  With  the  definitions 
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Table  1.4.  W  and  2/(W  +  1)  for  c/a  =  5,  o2>/avin>  =  1  x  106,  and  various  values 
of  Mq.  The  plasma  cloud  is  assumed  to  be  a  spheroid  with  symmetry 
transverse  to  the  magnetic  field. 


calculated.  The  results  are  shown  in  Table  1.4.  It  is  evident  from  the  table 
that  parallel  electron  currents  substantially  contribute  to  overall  quasi¬ 
neutrality  and  reduce  the  contribution  of  induced  Pedersen  currents  to  the  total 
perpendicular  current  in  the  ellipsoidal  plasma  cloud.  For  Mg  =  100,  cEox/B  =  1 
x  104  cm/s,  and  T.  (ion  temperature)  =  0.1  eV,  then  (1.55)  and  Table  1.4  give  kt 
=  0.2  km.  The  corresponding  wavelength  is  2ir/kj.  =  0.3  km.  This  wavelength  is 
consistent  with  the  perpendicular  scale  sizes  of  "frozen"  striations  [Linson  and 
Meltz,  1972]. 

Now  Linson  and  Meltz  [1972]  invoke  parallel  electron  currents  and 
coupling  to  conjugate  zones  in  order  to  account  for  the  observed  motion  of 
unstriated  barium-ion  clouds.  Effectively,  this  permits  Nc/Nq  (i.e.,  the  ratio 
of  height-integrated  cloud  to  ambient  electron  densities)  to  be  decreased 
relative  to  the  situation  when  parallel  electron  currents  and  coupling  to 
conjugate  zones  is  neglected.  Now  Nc/Ng  £  10  when  parallel  electron  currents 
and  coupling  are  considered  by  Linson  and  Meltz  [1972].  If  Nc/Ng  is  identified 
with  W  in  (1.53)  then  good  qualitative  agreement  is  obtained  between  Linson  and 
Meltz  [1972]  and  Table  1.4  (i.e.,  W  5.  10)  even  though  Linson  and  Meltz  [1972]  do 
not  explicitly  include  finite  parallel  conductivity  in  their  expressions  and 
this  present  paper  does  not  consider  an  inhomogeneous  background  ionosphere. 
Evidently,  parallel  electron  currents  play  an  essential  role  in  the  dynamics  of 
ionospheric  barium  clouds  by  substantially  reducing  the  magnitude  of  induced 
perpendicular  electric  fields  in  the  clouds. 


b. 


Temporal  Variations 


In  the  most  general  sense,  including  temporal  variations  of  external 
currents  and  induced  electric  fields,  the  induced  x-directed  current  within  the 
plasma  cloud  is  determined  from  (1.5b),  (1.23a),  and  (A4)  of  Appendix  A  to  be 
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ambient  ions.  Equation  (1.58)  differs  from  the  previously  hypothesized  view 
that  striation  "freezing"  indicates  the  balance  between  bifurcation  and  cross¬ 
field  diffusion  [McDonald,  et  al.,  1981;  Vesecky,  et  al.,  1980].  In  contrast  to 
the  earlier  hypothesis,  enhancements  in  the  cross-field  diffusion  need  not  be 
invoked  in  the  derivation  of  (1.58)  and  so  with  (1.58)  one  can  easily  account  for 
the  long-term  persistence  of  striations.  Moreover,  when  (1.58)  is  satisfied,  ion 
velocities  and  ion-viscous  forces  are  both  minimized  within  the  bulk  of  a 
striation. 

Figure  1.5  suggests  that  the  formation  of  striations  can  proceed  if 
plasma  condition  corresponds  to  points  above  the  curve.  This  is  illustrated  by 
trajectory  a  in  the  figure.  Conversely,  the  formation  of  striations  is  not 
likely  to  proceed  for  plasma  conditions  corresponding  to  points  below  the  curve 
in  Figure  1.5.  However,  for  these  plasma  conditions  it  still  can  be  hypothesized 
that  plasma  processes  like  ion-neutral  diffusion  may  tend  to  move  the  total 
plasma  condition  along  trajectory  b.  A  key  point  regarding  both  trajectories  a 
and  b  is  that  the  value  of  W,  attained  once  the  curve  is  reached,  must  be  larger 
than  the  initial  value  for  W.  This  is  consistent  with  the  idea  that  the 
formation  of  striations  must  be  a  manifestation  of  the  elongation  and  deformation 
of  a  plasma  cloud  in  response  to  an  applied  force.  Once  the  conditions  of  the 
plasma  cloud  are  such  that  a  point  on  the  curve  in  Figure  1.5  is  reached,  it  is 
assumed  that  the  plasma  condition  is  such  as  to  remain  on  the  curve.  Further 
bifurcation  is  not  possible  since  an  increase  in  W  is  inconsistent  with  a  smaller 
value  of  L,  as  indicated  by  trajectory  c  in  Figure  1.5.  Hence,  striations  are 
"frozen."  Ambipolar  diffusion,  ablation,  recombination,  and  etc.  decrease  W 
over  the  long  term  and  consequently  slowly  cause  the  plasma  state  to  go  along 
trajectory  d  in  Figure  1.5.  The  perpendicular  scale  size  of  striations  slowly 
decreases  along  this  trajectory.  From  the  standpoint  of  power  spectral 
densities,  the  outer  scale  length  (i.e.,  k^)  approaches  the  inner  scale  length 

[e.g.,  in  (D9)  of  Appendix  D]  [Vesecky  et  al.,  1980]. 

At  a  time  of  roughly  ten  minutes  after  release  at  approximately  200  km 

altitude  might  reasonably  be  expected  to  have  a  length  —  20  km  along  the  magnetic 

field  and  a  radius  2  km  transverse  to  the  magnetic  field  [Linson  and  Baxter, 
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1977].  Hence,  c/a  -  5.  If  oz>/a>v-n>  ~  1  x  10  then  ct/at  =  2.5  x  10  .  For  a 

spheroidal  plasma  cloud,  entry  4  of  Table  1.2  then  permits  W  and  2/ ( W  +  1)  to  be 
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The  curve  illustrating  the  marginal  condition  for  the  formation  of  striations  can 
be  determined  as  a  function  of  W  and  the  effective  Reynold's  number 


R  =  ~  .  (1.56) 
D 

In  (1.56)  the  following  symbols  are  used:  V  (effective  plasma  velocity),  L 
(effective  perpendicular  scale  length),  and  D  (effective  cross-field 
diffusivity) .  If  in  ( 1 . 54) - ( 1 . 56)  the  following  identifications  are  made 
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then 


W  +  1 

R  =  ,  (1.58) 

is  the  condition  for  marginal  stability.  This  expression,  which  is  plotted  in 
Figure  1.5,  does  not  have  an  explicit  dependence  on  Mq  but  is  nevertheless  valid 
only  for  Mq  >  1.  Equation  (1.58)  is  appropriate  to  lower  values  of  the  plasma 
density  for  which  ion  viscosity  is  unimportant. 

For  a  static  situation,  2EQX/(W  +  1)  in  (1.55)  can  be  interpreted  to  be 
the  ambipolar  electric  field  consistent  with  ion-neutral  diffusion  and  the  scale 
length  k~*.  This  suggests  a  physically  meaningful  interpretation  for  (1.58)  and 
the  curve  in  Figure  1.5.  The  equation  and  figure  indicate  the  balance  of  two 
effects:  bifurcation  and  the  attainment  of  thermodynamic  equilibrium  for  ions. 
Consider  an  ionospheric  barium  cloud.  By  thermodynamic  equilibrium,  we  mean  that 
(1.55)  is  satisfied  so  that  the  applied  electric-field  force  acting  on  cloud  ions 
is  exactly  balanced  by  the  pressure  forces  associated  with  the  cloud  ions  [cf, 
(E6)  in  Appendix  Ej.  Hence,  the  presence  of  the  cloud  does  not  result  in 
additional  ion  currents  within  the  frame  of  reference  where  the  effective 
electric  field,  2EQX/(w  +  1),  is  zero.  Stated  another  way,  the  current  carriers 
in  the  indicated  frame  of  reference  can  be  conceptualized  to  be  electrons  and 


strongly  coupled  [Perkins  et  al.,  1973]  and  undoubtedly  have  complicated  three- 
dimensional  dependences.  Indeed,  (1.21a-b)  demonstrate  that  induced  electric 
fields  can  have  strong  spatial  variations  for  even  the  relatively  simple  case  of 
a  uniform  ellipsoidal,  plasma  cloud  surrounded  by  a  uniform  plasma.  Furthermore, 
plasma  "equilibrium"  is  not  an  obvious  nor  proven  requirement  for  the  formation 
of  striations  L.  Wittwer,  private  communication,  1983  .  An  important 
unresolved  issue  is  the  understanding  of  how  three-dimensional  effects  impact 
the  basic  plasma  state  upon  which  striations  evolve. 

The  real  part  of  the  frequency  in  (1.54)  is  a  direct  consequence  of  the 
contribution  of  ion-neutral  diffusive  currents  to  overall  quasineutrality  and 
the  neglect  of  ambipolar  electric  fields  in  the  basic  plasma  model  upon  which  the 
striations  develop.  Now  if  ion-neutral  diffusive  currents  are  important. 
Appendix  E  demonstrates  that  the  contribution  of  electric  fields,  with 
associated  Pedersen  currents,  to  the  maintenance  of  quasineutrality  is 
correspondingly  reduced.  Consequently,  elongation  and  bulk  plasma  motion  in  the 
direction  of  the  force  is  also  inhibited.  Although  ion-neutral  diffusion  has  not 
been  considered  in  the  quantitative  evaluation  of  W  in  (1.52),  it  is  evident  from 
a  qualitative  standpoint  that  ion-neutral  diffusion  can  effectively  result  in 
larger  values  of  W  in  (1.53). 

It  is  assumed  that  bifurcation  can  proceed  only  if  Re(to)  <  Im(uo) 
[Sperling  and  Glassman,  1983;  Sperling,  1983a].  Presumably,  for  Re(w)  >  Im(u)) 
ion-neutral  diffusion  plays  a  dominant  role  in  maintaining  quasineutrality  and 
electric  fields  are  not  strongly  induced  in  the  nonlinear  striations.  Hence,  the 
tendency  of  striations  with  Re(u>)  >  Im(u)  to  elongate  a  plasma  cloud  and  cause 
plasma  motion  in  the  direction  of  the  applied  force  is  reduced. 

If  striations  can  only  develop  for  perpendicular  wave  numbers  such  that 
Re(ui)  5>  Im(u)),  the  permissible  perpendicular  wave  numbers  for  striations 
satisfy 
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For  lJx<totl  <<c  ! Jox< I  ^  1S  necessary  that  W  »  1  or  equivalently 
MQhx(Q)  »  1.  Clearly  Mq  »  1  is  not  a  sufficient  condition  for  the  total 
x-directed  current  in  the  ellipsoid  to  be  much  smaller  than  |  j'  . I.  This  is 
clearly  demonstrated  by  Figure  1.4  which  is  a  plot  of  W  versus  cf/a£  for  a 

22  1  21 2 

spheroid.  The  figure  shows  that  W  -►  1  for  c^/a^.  -*■  0,  but  W  -*•  Mq  for  ct/at  -►<*>. 
If  jx<  tQt  and  jQX<  are  ion-Pedersen  current  densities  then 
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In  (1.53)  the  total  x-directed  electric  field  in  the  ellipsoid  is  represented  by 
Ex<,tot‘ 

The  formation  of  striations  is  a  manifestation  of  the  tendency  of  plasma 
clouds  to  elongate  and  move  in  the  direction  of  an  applied  force.  It  follows 
that  with  the  development  of  striations  the  effective  overall  value  of  W  for  a 
plasma  cloud  must  increase.  Under  such  circumstances,  (1.51)  and  (1.53)  indicate 
that  jx<  tQt  and  Ex<  tot  are  appropriately  reduced. 

The  linear  frequency  of  collisional  E  x  B  or  Rayleigh-Taylor  instabil¬ 
ities  is  determined  by  combination  of  (1.21)  of  Sperling  [l983a]  and  (1.51)  or 
(1.53).  Specifically, 
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In  (1.54)  the  following  symbols  are  used:  to  (frequency),  ky  (wave  number 
transverse  to  the  ambient  magnetic  field  and  the  density  gradient),  (scale 
length  of  the  perpendicular  density  gradient),  and  P..  ( ion-gyroradius) .  The 
expression  assumes  vin>/fi-j>  =  vin</^i<  anc*  ^i>  =  To  account  for  other 
sources  of  current,  besides  electric  fields,  EQX  in  (1.54)  may  be  more  generally 
replaced  by  J0 x</^<- 

The  plasma  model  used  to  derive  (1.54)  includes  neither  ambipolar 
electric  fields  nor  parallel  electron  currents  as  both  these  phenomenon  can  be 
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assumed  to  be  uniform  in  the  perpendicular  direction.  Plasma  properties  are 
constant  for  |z|  <  zQ  and  |z|  >  zQ.  The  only  discontinuity  in  plasma  properties 
occurs  at  |z|  =  zQ.  For  simplicity,  it  is  assumed  here  that  the  plasma  geometry 
is  completely  symmetric  about  z  =  0. 

The  modes  are  assumed  to  have  the  plane  wave  form,  expfi^y  +  k>(<)2)  + 
Yt]  with  k^,  y,  Y ,  and  t  being  the  perpendicular  wave  number,  coordinate  in  the 
y-direction,  growth  rate,  and  time,  respectively.  Throughout  it  is  assumed  that 
the  subscripts,  >(<),  represent  quantities  appropriate  to  |z|  >(<)  Zq,  and  so 
k >(< )  represent  the  parallel  wave  number  outside  (inside)  the  unstable  plasma. 

In  general,  ionospheric  plasma  clouds  evolve  dynamically  as  a  result  of 
external  forces  and  diffusion  processes,  and  so  it  is  assumed  that  Y  »  1/t, 
where  t  is  the  time  scale  for  bulk  plasma  parameters  to  change. 


Sperling  [1983d]  derives  a  dispersion  relation  for  the  Rayleigh-Taylor 
and  E  x  B  instabilities  which  includes  the  contributions  of  parallel  currents, 
parallel  electric  fields,  and  electromagnetic  effects  in  a  medium  which  is 
uniform  along  the  magnetic  field.  This  dispersion  relation,  (15)  of  Sperling 
[l983d],  is 

(oj  -  u)q  +  ivi )  [l  -  TgUi )]  +  <^r0(Xi)  ^i 

u)  -  u)g  +  iv^  r  -  TQ(Xi  )J  w 
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The  following  symbols  are  also  used  in  writing  (2.1)  and  (2.2):  v.  (ion 
collision  frequency),  Iq  (modified  Bessel  function  of  the  first  kind  and  order 
zero),  Te  (electron  temperature),  (ion  temperature),  (Boltzmann  constant), 
mi  (ion  mass),  mg  (electron  mass),  e  (electron  charge),  ng  (electron  density), 
k |j  ( =  k>(<)>  parallel  wave  number),  u>  (=  iY,  frequency),  and  ve  (effective 
collision  frequency  for  electron  motion  along  the  geomagnetic  field).  Ions  are 
assumed  to  be  singly  ionized. 

The  parameter  Vq  is  the  zero-ordei ,  y-directed  velocity  of  individual 
ions.  For  example,  Vq  is  the  ion-Pedersen  drift  for  E  x  B  gradient  drift  modes 
[Ossakow  et  al.,  1979],  the  ion-polarization  drift  for  Rayleigh-Taylor  modes 
driven  by  ion  inertia  [Sperling,  1982b],  and  the  ion-gravity  drift  for  Rayleigh- 
Taylor  modes  driven  by  gravity  [Kadomtsev,  1965].  Expressions  for  Vq  are 
explicitly  written  and  described  in  the  references. 

We  now  consider  the  limit  of  collisional  modes  with  relatively  long 
wavelength,  i.e., 

v.  »  M  »  |u,t(,  |w*(,  |oi0| ,  (6 1 ;  Xi  «  1  •  (2.3) 
Equation  (2.1)  becomes 
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With  the  replacement  of  w  by  iY,  (2.4)  can  be  solved  for  6.  The  result  is 


the  maximum  growth  rate  for  the  pure  flute  instability.  Equation  (2.5)  implies 
forms  for  the  parallel  wave  numbers  in  the  interior  and  exterior  regions 
specified  in  Figure  2.1.  Specif ical ly. 
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Here  vA  is  the  Alfven  speed.  Because  of  the  use  of  (2.3),  (2.7b)  does  not  .depend 
upon  vQ  in  the  external  region  with  \z\  >  zQ. 

Equation  (13b)  of  Sperling  [1983d] ,  (2.3),  and  Faraday's  law  imply  that 
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In  (2.9)  modal  fields  are  represented  in  the  following  way:  B^x  (x-directed 
magnetic  field),  Elz  (z-directed  electric  field),  and  Ely  (y-directed  electric 
field) . 
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The  assumed  plasma  model  is  symmetric  about  z  =  0  and  (2.7b)  implies  that 
evanescence  away  from  the  plasma  cloud  is  the  appropriate  boundary  condition  for 
|z|  >  Zg.  The  following  explicit  forms  are  suggested  for  the  modal  electric  and 
magnetic  fields: 

Ely(y,z,t)  =  Ely<[exp(ik<z)  ±  exp(-ik<z)] 

H(z0  -  |z|)  +  Ely>  expf-lkju  -  Zg)] 

H(z  -  z0)  1  Ely>  exp[jk>| (z  +  z0)] 

H(-z  -  Zg)  expO^y  +  Yt)  (2.10a) 

Elz(y,z,t)  =  Elz<[exp(ik<z)  *  exp(-ik<z)] 

H ( Zg  -  I z | )  +  Elz>  exp  [-jkj  (z  -  Zg)] 

H(z  -  z0)  ?  Elz>  exp[|k>|  (z  +  zQ)] 

H(-z  -  Zg)  |  exp(ikxy  +  Yt)  (2.10b) 

Blx(y,z,t)  =  jBlx<[exp(ik<z)  +  exp(-ik<z)] 

H ( z g  -  |z| )  +  Blx>  exp [-J k>)  (z  -  z0)] 

H(z  -  Zg)  ?  Blx>  exp [|  k>|  (z  +  zQ)] 

H(-z  -  Zg)  |  exp(ikiy  +  Yt)  .  (2.10c) 

In  (2.10a-c)  the  Heaviside  function,  H,  is  equal  to  1  for  positive  argument  and  0 
for  negative  argument.  Consistent  with  (2.9a-b),  the  form  of  (2.10a-c)  indicates 
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that  Elz(y,z,t)  and  Blx(y,z,t)  are  odd  (even)  about  z  =  0  if  Ely(y,z,t)  is  even 
(odd)  about  z  =  0. 

To  complete  the  specification  of  the  eigenmodes,  it  is  necessary  to 
invoke  boundary  conditions  at  |z|  =  Zq.  In  particular,  tangential  electric  and 
magnetic  fields  are  continuous  across  a  discontinuity  in  dielectric  properties 
[Jackson,  1962;  Allis  et  al.,  1963;  Paoloni,  1975;  Sperling,  1983c].  For  the 
eigenmodes  with  Ely(y,z,t)  even  about  z  =  0,  it  follows  from  (2.9)-(2.10)  that 

Ely>  =  2Ely<  cos(k<20)  (2-Ua) 

“87  Ely>  =  2  £  Ely<  sin(k<zo)  *  (2-llb) 


Hence, 

^Zq  =  tan"1(|k>l8</k<8>)  +  pfT  (2.12) 

« 

with  p  a  nonnegative  integer  for  positive  k<.  Similarly  for  the  eigenmodes  with 
Ely(y,z»t)  odd  about  z  =  0, 

Ely>  '  21Ely<  sin(k<zQ)  (2.13a) 

-  Ely>  *  21  Eiy<  C0s(k<zo>  >  (2'13t>) 

and  so 

^Zq  =  -  tan’^k,  8^  |k>|  B<)  +  pn  (2.14) 


with  p  a  positive  integer  for  positive  k<.  For  a  specified  value  of  Zq,  the 
parallel  wave  numbers,  growth  rates,  and  modal  fields  are  completely  specified  by 
(2.7),  (2.9),  and  (2.12)  or  (2.14). 


Now  tan~^  is  bounded  by  0  and  n/2  and  so  the  lowest  order  eigenmode  is 
determined  by  (2.12)  with  p  =  0.  In  particular. 


Mq  =  tan”1  ( |  kj  8</k<6>) 


(2.15) 


If  | k> B</k< 6> |  «  1  then  it  follows  that  |  ^Zq)  is  also  much  smaller  than  1. 
Equation  (2.7)  indicates  that  the  inequality  has  a  tendency  to  be  better 
satisfied  if  the  density  of  an  unstable  plasma  cloud  is  much  larger  than  the 
density  of  the  stable  background.  The  examples  presented  in  the  subsequent 
section  satisfy  the  inequality.  Equations  (2.7),  (2.8),  and  (2.15)  now  give 


<V6>>' 


(2.16) 


For  Y  real  it  is  evident  from  (2.16)  that  0  <  Y  <  Yq.  By  squaring  both  sides  of 
(2.16)  it  is  possible  to  obtain  a  cubic  equation  of  the  form 


0  =  bu3  +  (1  -  a  -  2b)  u2  +  (b  -  2)  u  +  1 


(2.17) 


u  =  Y/Yr 


2  2  9 

a  !!i>  i_  h  -  “e>  Y0 

ve>  me  N*  kf  ’ 

i<  c  1  l  e> 


(2.18) 


The  parameter  "a"  indicates  that  the  total  integrated  electron  density  in  the 
unstable  cloud  (i.e.,  Nc  s  2ne<ZQ)  is  a  significant  quantity  [S.  Zalesak, 
private  communication,  1983]. 

It  is  instructive  to  examine  (2.17)  in  several  limits.  For  b  -*  0, 
(2.7)-(2.9)  indicate  that  the  modal  electric  fields  are  becoming  increasingly 
electrostatic  while  the  modal  magnetic  fields  are  becoming  smaller.  Equation 
(2.17)  then  can  be  reduced  to 


\  •*-  >  - 
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r 


0  s  (1  -  a)u  -  2u  +  1 


(2.19) 


with  appropriate  solution 


1  -  a 


1  -  a 


The  asymptotic  forms  for  (2.20)  are 


u  «  a-1'2  , 


u  «  1  -  a1/2  , 


(2.20) 


for  a  »  1 


for  a  «  1 


(2.21a) 


(2.21b) 


Equations  (2.20)  and  (2.21)  have  the  physically  attractive  feature  that  u  -►  0 
for  very  short  plasma  clouds  but  approaches  1  for  long  plasma  clouds.  In  the 
electrostatic  limit  (2.7)  becomes 
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(2.22a) 


(2.22b) 


In  the  electrostatic  limit  both  the  growth  rates  and  the  parallel  wave  numbers 
depend  on  k  .  Note  that  (2.22b),  which  is  independent  of  Y,  implies  that 
eigenmodes  must  become  increasingly  localized  along  the  magnetic  field  as  k^ 
increases. 

In  the  inductive  limit  with  b  -►«,  (2.17)  can  be  rewritten  in  the 
approximate  form 


0  =  u2  -  (2  +  a/b)  u  +  1 


(2.23) 


with  appropriate  solution 


(2.24) 


Asymptotic  forms  for  this  expression  are 

u  %  b/a  ,  for  a/b  »  1  (2.25a) 

u  a  l  -  (a/b)*^  ,  for  a/b  «  1  .  (2.25b) 

Again  the  asymptotic  expressions  for  u  have  the  pleasing  form  of  being  very  small 
for  sufficiently  short  clouds  and  approaching  1  for  sufficiently  long  clouds.  In 
the  inductive  limit,  (2.7a-b)  become 


vi<Y 


(2.26a) 
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(2.26b) 


Growth  rates  and  eigenmodes  do  not  depend  on  in  the  inductive  limit.  Also, 
(2.9a-b)  indicate  that  consistent  with  j b |  -*•  ®,  | E^2 1  -*■  0  while  | B^x |  becomes 
independent  of  k^ . 

2.3  RESULTS 

The  plasma  model  outlined  in  section  2.2  is  rather  simple  and  admittedly 
does  not  include  all  the  intricacies  and  complicated  inhomogeneities  of  the  real 
ionosphere.  Nevertheless,  the  model  does  show  that  growth  rate  is  a  function  of 
the  field-line  electron  density  of  the  plasma  cloud  (i.e.,  Nc).  In  the  inductive 
limit,  growth  rates  and  parallel  wave  numbers  are  not  a  function  of  the 
perpendicular  wave  number.  In  the  electrostatic  limit,  growth  rates  and  parallel 
wave  numbers  are  strong  functions  of  the  perpendicular  wave  number. 


As  quantitative  examples,  consider  Figures  2.2a-b  which  are  plots  of  k 

1  12-2  1 
(km"A)  versus  Y/YQ  for  different  values  of  Nc  (10  cm  ).  Figures  2.2a-b  assume 

YQ  =  0.01  s"1  and  YQ  =  0.1  s"1,  respectively.  The  plots  assume  ng>  = 

1  x  105  cm"3,  Ti<  =  Tj  =  Te<  =  Tg>  =  0.1  eV,  and  B  =  0.4  G.  The  ions  for 

| z |  <  Zq  are  assumed  to  be  barium  (i.e.,  m^  =  137  m  ,  with  rrip  the  proton  mass) 

and  the  ions  for  |z|  >  Zg  are  assumed  to  be  air  with  m^  =  20  m  .  The  mass  of  an 

atmospheric  neutral  particle,  mo <:(>)*  a^so  assumec*  to  be  ^0  mp*  The  neutral 

density,  n0<  =  nQ>  =  8.3  x  109  cm-3,  is  appropriate  to  an  altitude  of  ~200  km 

[Knapp  and  Schwartz,  1975].  The  collision  frequencies  are  determined  from 

[Braginskii,  1965;  Kilb,  1977] 


i<  =  2*5  x  10  mgng 

(2.27a) 

^  =  2.5  x  1033  mgng 

(2.27b) 

e<( >)  =0*51  vei<(>)  +  1*9  x  10"8  n0 

(2,27c) 

:1<  =  °*51  vei<  +  1*9  x  10"8  n0 

(2.27d) 

ei<( >)  =  1,1  x  10  3  ne<( >) 

(2.27e) 

In  (2.27)  all  parameters  are  in  cgs  units. 

Figures  2.2a-b  show  the  following  general  tendencies.  The  growth  rate, 
Y,  becomes  larger  as  the  perpendicular  wave  number,  k  ,  increases  and  as  Nc 
increases.  The  modes  become  more  inductive  with  I  3Y/3k  |  -*■  0  as  Nc  increases 
and  k  becomes  smaller.  The  inductive  nature  of  modes  is  enhanced  if  the  maximum 
possible  growth  rate  for  instability,  YQ,  increases.  For  larger  values  of  k± 
(i.e.,  ki  -*•  10  km“^),  modes  become  highly  electrostatic.  In  accordance  with 
(2.20),  Y/Yg  then  is  independent  of  Yq  and  so  the  corresponding  curves  in 
Figures  2.2a-b  become  essentially  the  same.  Equation  (2.16)  is  a  good  approxi¬ 
mation  to  (2.15)  for  Y/Yg  <  0.95. 
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Consider  a  plasma  which  is  simultaneously  steepening  perpendicular  to 
the  magnetic  field  and  elongating  along  the  magnetic  field.  If  ionization, 
recombination,  diffusion,  and  convection  can  be  neglected,  a  reasonable 
assumption  is  that  Nc  is  a  constant  as  the  unstable  plasma  cloud  elongates 

[S.  Zalesak,  private  communication,  1983].  Earlier  in  the  steepening,  Yn  might 

-1  u 
be  0.01  s  and  Figure  2.2a  is  appropriate.  The  figure  clearly  shows  the 

tendency  for  growth  rate  to  be  a  strong  function  of  ,  particularly  for  smaller 

values  of  N  .  Parallel  currents  and  parallel  electric  fields  tend  to  most 

strongly  short-circuit  perpendicular  electric  fields  associated  with  the 

unstable  modes  when  k  is  small.  It  is  evident  that  preferential  growth  occurs 

1 

at  larger  wave  numbers.  At  later  times,  when  Yg  might  be  ~0.1  s  ,  the  same 
basic  tendencies  still  remain,  but  growth  rates  are  a  somewhat  weaker  function  of 
k^  since  inductive  effects  are  more  important. 

When  E  x  B  gradient  drift  instabilities  are  the  primary  cause  of 

striations  in  barium  clouds,  the  onset  of  striations  in  barium  clouds  occurs' 

several  minutes  after  the  detonation  of  the  barium  release.  Specifically  for 

larger  barium  clouds,  like  Spruce  [Linson  and  Meltz,  1972],  the  onset  time  for 

striations  is  estimated  to  be  generally  greater  than  approximately  5  minutes 

[McDonald  et  al.,  1980],  Now  Linson  and  Meltz  [1972]  have  estimated  that- the 

Spruce  and  Olive  ion  clouds  had  an  approximately  30-40  km  extent  along  the 

magnetic  field  at  approximately  5  minutes.  This  parallel  extent  corresponds  to 

zn  15-20  km  for  Figures  2.2a-b.  If  n  =  1  x  107  cnf^,  then  N  =  3-4 

x  101,3  cm  .  For  the  growth  of  striations  to  be  significant,  it  is  assumed  that 

the  growth  rate  for  striations  must  be  greater  than  the  rate  at  which  the  bulk 

parameters  of  the  cloud  change,  but  finite  Nc  implies  a  growth  rate  which  is 

smaller  than  Yg.  Now  if  the  typical  time  scale  for  bulk  plasma  parameters  to 

change  is  t  =  300  s,  then  Y  »  1/t  implies  Y/Yq  »  0.3  and  0.03  in 

Figures  2.2a  and  2.2b,  respectively.  Figure  2.2a  indeed  shows  that  clouds  with 
13  -2 

Nc  £  2  x  10  cm  have  a  substantial  range  of  finite  perpendicular  wave 

numbers  with  sufficiently  rapid  growth  rates.  Figure  2.2b,  which  is  appropriate 

to  the  case  of  more  substantial  steepening  than  Figure  2.2a,  indicates  that 

13  -2 

Y  »  1/300  s  is  easily  satisfied  for  Nc  >  2  x  10  cm  .  Examination  of 

Figures  2.2a  and  2.2b  suggests,  that  in  evolving  plasma  clouds,  preferential 
striation  growth  should  occur  for  k^  finite  rather  than  for  k^  =  0. 


2.4  CONCLUDING  REMARKS 

A  simple  plasma  geometry  is  used  to  show  that  the  finite  field-line 
integrated  electron  density  of  ionospheric  plasma  clouds  plays  a  role  in  the 
delayed  formation  of  striations  by  tending  to  reduce  the  growth  rate  of  striation 
instabilities.  Taking  into  account  the  field-line  integrated  electron  density, 
the  onset  of  appreciable  striations  is  assumed  to  occur  when  the  growth  rate  of 
the  striation  instabilities  is  much  faster  than  the  rate  of  change  of  the  bulk 
plasma.  This  point  is  quantified  by  examination  of  Figures  2.2a  and  2.2b.  The 
finite  parallel  length  of  plasma  clouds  tends  to  favor  the  growth  of  striations 
with  finite  rather  than  with  zero  perpendicular  wave  numbers. 

The  results  presented  in  this  chapter  are  not  qualitatively  changed  by 
considering  more  complicated  and  realistic  plasma  geometries  than  the  one  used 
here.  However,  it  should  be  noted  that  warm  plasma  effects  (e.g.,  ambipolar 
diffusion),  not  considered  in  the  present  chapter,  tend  to  reduce  growth  rates 
for  sufficiently  large  k  [McDonald  et  al.,  198l] .  Indeed,  the  combination  of 
parallel  inhomogeneities  and  warm  plasma  effects  generally  causes  a  maximizing 
of  the  growth  rate  for  k  in  the  range  of  1  ~  10  km-1  [j.  L.  Sperling  and 
A.  J.  Glassman,  1983]. 
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3.  LATE-TIME  STRIATION  EIGENMODES  ALONG  THE  GEOMAGNETIC  FIELD 
3.1  INTRODUCTION 

Following  a  high-altitude  nuclear  burst,  plasma  is  injected  at  high 
speed  along  the  geomagnetic  field  into  the  upper  ionosphere  and  the  magneto¬ 
sphere.  The  resulting  nuclear  plumes  can  have  densities  orders  of  magnitude 
larger  than  under  normal  ambient  conditions.  In  fact,  in  view  of  their  large 
overall  spatial  extent  and  high  density,  nuclear  plumes  represent  a  markedly 
different  medium  than  spread-F  plasmas  or  barium  clouds  for  the  formation  of 
striations.  For  example,  it  has  been  reported  that  collisional  striation 
instabilities,  associated  with  ion-neutral  collisions,  can  be  expected  to  have  a 
much  shorter  extent  along  the  geomagnetic  field  in  the  high-altitude  nuclear 
environment  than  in  the  barium-cloud  environment  [Sperling,  1983e]. 

The  large  parallel  extent  and  high  density  of  nuclear  plumes  offers  the 
possibility  that  collisionless  fluid  instabilities,  driven  at  high  altitude,  may 
play  a  very  substantial  role  in  structuring  at  later  times  [Kilb,  1977].  Of 
course,  because  nuclear  plumes  extend  well  into  the  ionosphere,  the  possibility 
remains  that  collisional  instabilities,  associated  with  ion-neutral  collision, 
can  also  still  provide  .impetus  for  structuring.  By  contrast,  the  smaller  mass 
and  parallel  extent  of  barium  clouds,  likely  limits  their  structuring  to  solely 
collisionless  or  collisional  processes  depending  on  the  altitude  of  the  barium 
release.  For  barium  releases  at  an  altitude  of  approximately  200  km,  structuring 
is  probably  a  result  of  collisional  instabilities. 

In  nuclear  plumes,  the  curvature  of  the  geomagnetic  field  and  the  high 
velocity  (e.g.,  -  5  km/s)  at  which  plasma  is  ejected  along  the  geomagnetic  field 
suggests  that  centrifugal  forces  might  contribute  to  the  "effective  gravity," 
necessary  for  structuring  by  a  collisionless  Rayleigh-Taylor  mechanism  [Kilb, 
1977].  Assuming  a  dipolar  model  for  the  geomagnetic  field,  an  appropriate 
expression  for  gc  is  [Kilb,  1977]: 


A 


q  =  fl!iL  [l  +  cos2(9)]sin(8)  A  Bo 
c  R  [l  ♦  3  coS2(e)]3«  *  Bo 


(3.i)  ■: 
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(3.21).  In  particular,  4>(z)  is  required  to  decay  far  from  the  source  of  free 
energy  responsible  for  instability,  i.e., 

I  'l'  ( z-^±oo)|  -*■  0  .  (3.30) 

The  appropriateness  of  this  boundary  condition  is  shown  in  great  detail  in 
Sperling  [l983e]  to  mainly  result  from  the  tendency  of  striation  fields  to  become 
more  electrostatic  at  lower  altitudes.  Indeed,  the  numerical  results  presented 
in  Section  3.4  demonstrate  the  correctness  of  (3.30). 

3.3  APPLICATIONS  TO  MODELING 

The  ability  to  bifurcate  requires  that  the  relative  phases  of  parameters 

(e.g.,  density,  electric  field,  etc.)  in  the  parent  structures  be  identical  to 

the  relative  phases  in  the  daughter  structures  [Vesecky  et  al.,  1980] .  However, 

for  k^  >  k^  our  numerical  results,  using  the  assumed  plasma  models,  indicate 

that  the  modes  become  highly  dispersive  and  the  relative  phases  of  parameters 

become  sensitive  functions  of  the  perpendicular  wave  number.  Here  kt  is  the 

perpendicular  wave  number  at  which  the  real  part  of  the  frequency  is  equal  to  the 

growth  rate.  Moreover,  since  the  eigenmodes  become  more  localized  along  the 

«  # 
geomagnetic  field  as  perpendicular  wave  number  increases  and  concomitantly, 

eigenvalues  become  insensitive  to  parallel  variations  of  the  plasma,  the  density 

of  eigenmodes  and  eigenvalues  in  phase  space  increases  along  with  the  real 

frequency.  This  suggests  that  striations  are  limited  to  scale  sizes  larger  than 

the  wavelength  for  which  the  real  frequency  is  equal  to  the  growth  rate. 

The  coefficient  for  ion  diffusion,  resulting  from  ion-neutral 
collisions,  can  play  an  important  role  in  determining  k^,  especially  when 
unstable  plasmas  are  not  too  dense.  This  is  demonstrated  in  Appendix  G  which 
determines  k^  in  various  important  limits.  The  ion-neutral  diffusion  does  not 
directly  stabilize  the  instability  because  electrons  tend  to  diffuse  across  the 
magnetic  field  at  a  much  slower  rate  than  ions,  but  it  does  cause  the  instability 
to  acquire  a  significant  real  frequency  for  the  assumed  plasma  model. 

Although  the  real  part  of  the  frequency  may  be  sensitive  to  the  assumed 
theoretical  plasma  models  [S.  Zalesak,  private  communication,  1984],  the 
perpendicular  wave  number  determined  by  equating  the  real  and  imaginary  parts  of 
the  frequency  can  be  interpreted  in  a  way  which  is  independent  of  the  models.  In 
particular,  at  least  for  the  collisional  gradient-drift  instability,  it  can  be 


For  uj  =ss  iY  with  Y  real  and  Y  »  |co  J,  the  effective  resistivity  is  increased  by 
2  2  2  ® 

the  ratio  1  +  (oigY/k^c  v  )  from  the  purely  classical  electrostatic  value.  Since 


u  -  U) 


1  «  P  = 


(3.26) 


implies  an  important  contribution  from  the  electromagnetic  correction,  it  is 

2  2  * 

evident  that  for  fixed  k.  ,  c  ,  u,  and  w  the  electromagnetic  correction  has  the 

i  e  2 

greatest  impact  at  higher  altitudes  where  <d  fa  is  maximum.  At  lower  altitudes 

2  c  G 

uig/v  is  reduced  by  electron-neutral  collisions  [Goldman  et  al.,  1976]  and  the 
electromagnetic  correction  is  less  significant.  It  is  interesting  to  note  that 
turbulence  and  double  layers,  which  increase  the  effective  value  of  v  ,  would 
also  reduce  the  importance  of  electromagnetic  effects  [Block,  1972; 
Papadopoulos,  1977;  Sperling,  1982b]. 

Consider  the  neglect  of  ion-viscosity,  ambipolar-diffusion,  and  ion- 
polarization  corrections  in  (3.21)  and  (3.22).  If  in  addition. 


1  » 


kia<BTeo 


(3.27) 


and  asymptotically  electrostatic  modes  are  assumed,  i.e.. 


ue  w  '  ^e 


(3.28) 


then  (3.21)  and  (3.22)  reduce  to  (13)  of  Goldman  et  al.  [1976].  Equation  (3.27) 
is  generally  easily  satisfied  for: 


Vi  I  •  l“ll«  M<VTeo> 


(3.29) 


To  completely  determine  and  the  other  associated  striation  parameters, 
it  is  necessary  to  specify  the  boundary  conditions  imposed  on  the  solution 
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In  (3.21) 


u>  2 

i  ~2~2  '  we)  " 


,  B(z)  = 


dA(z) 

dz 


(Vc  - 

’eo  “Veo  1 


(3.22) 


Equation  (3.21)  is  the  dispersion  relation  solved  in  this  work. 

With  the  form  of  4*  determined  from  (3.21),  n^,  n.j,  B^,  E^,  and  vg^z 
are  determined  from  (3.16),  (3.17a),  and  (3.17c-e).  In  addition,  the  parallel 
current  associated  with  the  modes,  Je^z»  is  determined  from: 


'■’elz  ~  eneovelz 


(3.23) 


If  ion-ion  diffusion  is  neglected,  (3.21)  can  be  easily  generalized  to 
multiple  ion  species  by  redefining  a  such  that: 


a  =  . 
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(3.24) 


with  the  summation  over  ion  species  and 


(3.25) 


2  *  2  2 

The  electromagnetic  correction  to  (3.21)  is  localized  to  u)g(co  -  w^/k^c 
in  the  denominator  of  A(z)  in  (3.22).  From  the  mathematical  standpoint,  the 
electromagnetic  correction  acts  to  effectively  modify  the  plasma  resistivity. 


76 


In  (3.17a-f) 
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The  second-order  differential  equation  can  be  written  in  a  somewhat 
simpler  form  with  the  definition: 
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Hence, 
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and 


0  =  |a(z)  +  B(z)  jj-  +  C(z)l  * 


(3.21) 
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In  writing  (3.15d),  the  displacement  current  has  been  neglected.  This  assumption 
is  equivalent  to  quasi-neutrality,  i.e.. 


nel  =  nil  • 


(3.16) 


The  combination  of  (3.12)  and  (3.15a-d)  gives: 


nel  =  aneoEly 


(3.17a) 
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velocity  vgj2  in  response  to  the  modes,  is  determined  from  the  linearized 
momentum  equation: 


0  <BTeoneo  dz  (ne)  "  eneoElz  '  meveneovelz 


+  eneoveoy®lx‘/c 


(3.12) 


with  [Braginskii,  1965] 


ve  =  ven  +  0,51  vei  *  veoy 


c  kbt  dna„  , 

B  eo  eo  1 

eB0  dx  neo 


(3.13) 


Here  veoy  is  the  zero  order  electron  diamagnetic  drift  and  B^x  is  the  x-directed, 
modal,  magnetic  field.  Hence, 
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Next  we  consider  the  ion  and  electron  continuity  equations,  Faraday's 
law,  and  Ampere's  law: 

0  -  -iunu  ♦  nilui0y  +  h  <n1ov11x>  +  iki"lov11y  >  (3'15a> 
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In  (3.9),  w.  is  the  ion-diamagnetic  drift  frequency.  It  follows  that: 
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In  (3.10a-b): 


(3.10c) 


(3. lOd ) 


1  +  v2/q2 


(3.11) 


The  second  term  in  (3.10b)  is  an  approximation  which  is  reasonable  provided 
ven  >:>  vei  when  a  >:>  Under  such  circumstances,  the  contribution  from  dif¬ 
fusion  due  to  electron-neutral  and  ion-neutral  collisions  is  much  larger  than  the 
diffusion  resulting  from  electron-ion  collisions. 

In  addition  to  the  perpendicular  electron  currents,  parallel  electron 
currents  contribute  to  the  growth  of  modes.  Parallel  ion  currents  are  much  less 
important  and  are  neglected  here  [Sperling,  1983d].  The  parallel  electron 


,•»  * *  -  ”  *"  .  *  _  •_  •  -  •* 
V  V  v.  *  •*/-* 


72 


«,  •  .  *«  •  7717  W  3?  • 


a1  neo 


3.5  x  1CT  T 


5  t3/2  » 


a'  =  23.4  -  1.15  1og10(neQ)  +  3.45  log1Q  (Te<J)  , 


%  8  x  10'8  Te‘64  "no  • 


v.  =  2.5  x  10  m  n 
in  n  no 


for  air  ions 


(3.7) 


All  parameters  are  in  cgs  units  except  for  temperatures  which  are  in  units  of  eV. 

-12 

The  Boltzmann  constant  is  <g  **  1.6  x  10  ergs/eV.  In  (3.7),  A,  mn,  and  nnQ 
represent  the  ion  atomic  mass,  neutral  particle  mass,  and  the  neutral  number 
density,  respectively.  The  expression  for  \>^  includes  contributions  from  both 
the  ion-neutral  collisions  and  ion-ion  viscosity.  In  sufficiently  dense 
ionospheric  clouds,  the  ion-ion  collision  frequency,  v^,  can  be  larger  than  the 
i  on-gyrof requency . 

In  (3.6a-b), 


wo  "  kiuioy  *  uioy  ~  v.  B 


in  o  1  +  r). 


*  ^in  ”  ft.. 
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and  u-  is  the  y-directed  drift  velocity  of  individual  ions.  The  y-component  of 
oy  a  a 

the  effective  gravity  is  written  as  g  .  Here,  z  and  y  are  the  unit  vectors  in  the 

z-  and  y-directions,  respectively. 

Equations  (3.6a-d)  are  solved  with  the  assumption  that: 

ai  »  »  U0I  .  |w*l  »  Ven/mi  ’  mevei/mi 
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(3.9) 


The  linear  ion  and  electron  velocities  perpendicular  to  the  background 
magnetic  field  and  in  response  to  the  modes  can  be  determined  from  the  fluid 
equations: 


-  <»0)v11y  =  -VloVil  *  enioEly  '  eniovilxVc 
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In  (3.6a-d)  the  subscript,  o,  indicates  background  parameters  while  the  sub¬ 
script,  1,  indicates  the  response  to  the  modes.  Ion  quantities  are  represented 
by  the  subscript,  i,  and  the  subscript,  e,  is  used  for  electron  quantities.  The 
x-  (y-)  directed  components  of  vectors  are  represented  by  the  subscript  x(y). 
The  following  symbols  are  also  used  in  (3.6a-d):  <B  (Boltzmann  constant),  T 
(temperature),  n  (density),  E  (electric  field),  v  (velocity),  vi  (effective  ion 
collision  frequency),  vg1-  (electron-ion  collision  frequency),  and  vgn  (electron- 
neutral  collision  frequency).  The  magnitude  of  the  electron  charge  is  repre¬ 
sented  by  the  symbol,  e.  In  (3.6a-d)  explicit  forms  for  and  vgi  can  be  written 
in  the  following  way  [Braginskii,  1965;  Ki lb,  1977]: 
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DISPERSION  RELATION 


To  derive  the  linear  dispersion  relation  for  E  x  B  gradient  drift  and 
Rayleigh-Taylor  eigenmodes  along  the  geomagnetic  field,  it  is  sufficient  to 
consider  a  low  beta,  slab-geometry  _  plasma  model  with  a  z-directed  constant 
magnetic  field  with  strength,  BQ,  a  density  profile  which  depends  on  the 
x-coordinate,  and  a  y-directed  constant  ambient  electric  field  with  strength, 
E0y,in  the  neutral-wind  frame  of  reference.  In  the  neutral-wind  frame  of 
reference,  electrons  and  ions  have  an  x-directed  E  x  B  velocity  of  size  cE0y/BQ. 
Plasma  parameters  are  permitted  to  vary  in  the  z-direction.  The  collective 
oscillations  are  assumed  to  have  the  plane  wave  form,  exp[i(k^y  -ut)]  ,  with  k  , 
u>,  and  t  being  the  perpendicular  wave  number,  frequency,  and  time,  respectively. 
The  z-dependence  of  the  oscillations  is  determined  self-consistently. 

The  plasma  model  includes  neither  ambipolar  electric  fields  nor  parallel 
electron  currents  as  both  these  phenomena  can  be  strongly  coupled  [Perkins 
et  al.,  1973]  and  undoubtedly  have  complicated  three-dimensional  dependences. 
For  example,  parallel  variations  in  the  perpendicular  density  gradient  would 
imply  that  the  perpendicular  electrostatic  potential,  driven  by  ion-neutral 
diffusion,  must  be  a  function  of  the  coordinate  along  the  geomagnetic  field  and 
so  parallel  electric  fields  are  present.  Indeed,  (1.21a-b)  of  Chapter  1 
demonstrate  that  induced  electric  fields  have  strong  spatial  variations  for  even 
the  relatively  simple  case  of  a  uniform  ellipsoidal,  plasma  cloud  surrounded  by  a 
uniform  plasma.  Furthermore,  plasma  "equilibrium"  is  not  an  obvious  nor  proven 
requirement  for  the  formation  of  striations  [l.  Wittwer,  private  communication, 
1983].  An  important  unresolved  issue  is  the  understanding  of  how  three- 
dimensional  effects  impact  the  basic  plasma  state  upon  which  striations  evolve. 

Fluid  modes  with  frequencies  much  smaller  than  the  ion-gyrofrequency  are 
considered.  As  usual  with  a  fluid  treatment,  it  is  assumed  that  electron  and  ion 
mean  free  paths  for  collision  are  much  smaller  than  the  scale  lengths  for 
parallel  spatial  variations  and  electron-  and  ion-gyroradi i  are  much  smaller 
than  the  scale  lengths  for  perpendicular  spatial  variations  [Boyd  and  Sanderson, 


ion-gravitational  drifts  are  stabilizing.  Under  such  a  circumstance,  the 
possibility  for  instability  is  confined  to  the  ionosphere  where  the  ion-Pedersen 
drift  is  finite.  As  indicated  by  entry  4  in  Table  3.1,  there  is  stability 
everywhere  if  both  the  ion-Pedersen  drift  and  the  combination  of  ion-centrifugal 
and  ion-gravity  drifts  are  stabilizing. 

To  determine  the  relative  importance  of  collisionless  and  collisional 
instabilities  in  late-time  phenomenology,  a  second-order  linear  differential 
equation  is  derived.  The  equation  determines  eigenvalues  ( i . e . ,  temporal 
frequencies  in  complex  space)  and  the  associated  eigenmode  structure  of  the 
instabilities  along  geomagnetic  field.  The  differential  equation  is  solved  for 
parameters  and  geometry  appropriate  to  nuclear  plumes. 

It  is  demonstrated  that  for  a  specified  late-time  plasma  configuration, 
striations  tend  to  develop  with  larger  scale  sizes,  perpendicular  to  the 
geomagnetic  field,  at  higher  altitudes  than  at  lower  altitudes.  Electromagnetic 
effects  are  shown  to  be  important  for  modes  with  smaller  perpendicular  wave 
numbers  (i.e.,  5,  1  km"^)  and  ion-viscosity  is  important  for  modes  with  larger 
perpendicular  wave  numbers  (i.e.,  &4  km"^").  Analytic  criteria  are  given  for  the 
perpendicular  wave  numbers  of  the  instability  at  which  the  real  part  of  the 
frequency  is  equal  to  the  growth  rate  and  at  which  there  is  complete  stabiliza¬ 
tion.  It  is  argued  that  the  equating  of  real  part  of  the  frequency  and  the  growth 
rate  gives  a  meaningful  estimate  of  the  perpendicular  wave  number  at  which 
striation  "freezing"  occurs.  An  appropriate  algorithm,  suitable  for  implemen¬ 
tation  in  modeling  codes,  is  given  for  the  perpendicular  wave  number  at 
"freezing."  Simple  empirical  estimates  for  the  parallel  spatial  extent  of  the 
striation  density  fluctuations  are  described. 

This  chapter  is  divided  into  five  sections.  Section  3.2  derives  the 
differential  equation  which  is  the  dispersion  relation  for  the  striation 
instabilities  along  the  geomagnetic  field.  Section  3.3  gives  applications  for 
modeling  the  formation  of  striations  in  late-time  HANE  plasmas.  Section  3.4 
presents  results  appropriate  to  late-time  HANE  phenomenology.  Section  3.5  is  a 
summary  and  gives  concluding  remarks. 
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(3.4) 


In  (3.4)  is  the  scale  length  of  the  perpendicular  density  gradient.  At  lower 
altitudes  within  the  ionosphere  where  the  ion-neutral  collision  frequency  is 
large  compared  to  the  growth  rate,  the  characteristic  growth  rate  of  the 
instability  has  the  form  appropriate  to  collisional  gradient-drift  or 
collisional  Rayleigh-Taylor  modes  [Volk  and  Haerendel,  1971;  Perkins  et  al., 
1973]: 
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For  the  present  purpose,  it  is  sufficient  to  note  that  Y  «  '  while  Y  «  ^  . 

This  suggests  that  the  collisionless  instability  is  favored  by  long  perpen¬ 
dicular  scale  lengths  in  the  density  profile  (i.e.,  large  L^)  while  the  col¬ 

lisional  mode  is  favored  by  short  scale  lengths  in  the  density  profile  (i.e., 
small  L^).  Now  in  general,  the  actual  relative  significance  of  collisionless  and 
collisional  modes  depends  on  values  for  parameters  and  geometry.  For  example, 
consider  Ig^l  =  5  x  102  cm/s2  in  (3.4)  and  f9el/vin  *  1  x  104  cm/s  in  (3.5).  If 
Lx  =  2  x  106  cm  then  1.6  x  10"2  s_i  =  Y  >  YQ  *  5  x  10~3  s"1,  but  if 

=  1  x  10^  cm  then  7.1  x  10  2  s”^  =  Y  <  Yg  =  1  x  10”2  s-^.  Now  consider  |  g"Cgl 

=  1  x  10^  cm/s2  in  (3.4)  and  |g_|/v.  =  1  x  10^  cm/s.  It  then  follows  that 

_  c  g  i  n 

Y  >  Yq  for  Li  >  1  x  10°  cm.  These  simple  quantitative  examples  illustrate  the 
potentially  important  role  for  collisionless  instabilities  in  whatever 
structuring  occurs  in  late-time  nuclear  plumes. 

Examination  of  (3.2)  indicates  various  possibilities  for  instability 
depending  upon  the  sign  of  9^Cg*Vneo  and  ^’^eo’  These  possibilities  are 

summarized  in  Table  3.1.  The  emphasis  in  this  report  is  placed  on  the  case  when 

the  ion-Pedersen  drift  and  the  combination  of  centrifugal  acceleration  and 
gravity  permit  the  possibility  for  instability  at  both  high  and  low  altitudes. 
This  case  is  specified  as  entry  1  in  Table  3.1.  Entry  2  of  Table  3.1  illustrates 
that  instability  tends  to  be  confined  to  higher  altitudes  if  the  ion-Pedersen 
drift  is  stabilizing.  Conversely,  entry  3  in  Table  3.1  shows  the  case  when  the 
ion-Pedersen  drift  is  destabilizing  but  the  combination  of  ion-centrifugal  and 
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In  (3.1)  the  following  symbols  are  used:  R  (distance  from  the  center  of  the 
earth),  wM  (plasma  speed  along  the  geomagnetic  field),  0  (geomagnetic  colati- 

A  -► 

tude),  e^  (unit  vector  in  the  geomagnetic  eastward  direction),  BQ  (geomagnetic 
field) ,  and  BQ  (magnitude  of  the  geomagnetic  field). 

Away  from  the  geomagnetic  poles,  the  geomagnetic  field  has  a  component  in 
the  direction  of  the  geomagnetic  colatitude.  Hence,  the  component  of  gravity 
perpendicular  to  the  geomagnetic  field  can  contribute  to  structuring  by  a  col¬ 
lisionless  Rayleigh-Taylor  mechanism  [Kadomtsev,  1965;  Hudson  and  Kennel,  1975; 
Ki lb,  1977]. 

The  "effective  gravity,"  g"e,  which  drives  striation  growth  at  late  times 
is  permitted  to  consist  of  centrifugal,  gravitational,  and  neutral  wind 
contributions  and  so  g*e  has  the  following  form. 


9e  9cg  +  9E 


(3.2) 


9cg  =  9c  +  9i  • 


j  Bo  x  CEo 

in  Bo  Bo 


(3.3) 


In  (3  2)  and  (3.3)  the  symbols  are:  (ion-neutral  collision  frequency),  EQ 
(applied  electric  field),  c  (speed  of  light),  and  g^  (gravitational  acceleration 
perpendicular  to  the  geomagnetic  field).  Equation  (3.3)  is  appropriate  to  a 
frame  of  reference  in  which  a  uniform  neutral  wind  is  stationary.  Define  the 
electron  density  to  be  neQ.  Since  Rayleigh-Taylor  instabilities  require 
9cg’Vneo  <  9ravity  ten^s  to  destabilize  the  underside  of  nuclear  plumes  on 
lower  L  shells,  while  the  centrifugal  force,  resulting  from  the  curvature  in  the 
geomagnetic  field,  tends  to  promote  stabilization.  Conversely,  on  the  upper  side 
of  a  nuclear  plume,  on  higher  L  shells,  the  centrifugal  force  is  destabilizing 
while  gravity  is  stabilizing. 

At  high  altitudes  where  the  ion-neutral  collision  frequency  is  'small 
compared  to  the  growth  rate,  the  characteristic  growth  rate  of  the  instability 
has  the  form  appropriate  to  collisionless  Rayleigh-Taylor  modes  [Kadomtsev, 
1965:  Hudson  and  Kennel.  1975]: 


demonstrated  that  the  form  for  can  be  obtained  by  balancing  the  applied 
electric  field  (in  the  neutral -wind  frame  of  reference)  and  the  ion  pressure 
gradient  within  a  striation  [see  (G12)  in  Appendix  Gj.  Under  such  a 
circumstance,  it  is  not  necessary  for  there  to  be  induced  electric  fields  in 
order  for  current-density  fields  to  be  nondivergent.  In  fact,  the  form  for  kt 
suggests  striation  "freezing"  represents  the  balance  between  two  effects: 
(1)  the  tendency  for  bifurcation  and  (2)  the  tendency  of  ions  to  approach  a  state 
of  thermodynamic  equilibrium  with  applied  electric  fields.  When  ions  are  in 
thermodynamic  equilibrium  with  applied  electric  fields,  electrons  are  the 
primary  current  carriers  in  the  frame  of  reference  where  the  basic  plasma  state 
is  stationary.  Moreover,  ion  viscosity  and  ion  thermal  conductivity  drive  ions 
toward  thermodynamic  equilibrium  and  have  no  effect  on  ions  when  thermodynamic 
equilibrium  is  attained.  Because  the  form  and  magnitude  for  kt  appears  to  be 
reasonable  for  barium  clouds,  it  is  assumed  here  that  similar  calculations  for 
the  high-altitude  nuclear  environment  are  also  reasonable. 

If  4>(z)  can  be  approximated  as  having  a  plane-wave  dependence  on  z,  i.e., 
exp(ik^z),  then  by  the  replacement 

5?  *  (3-31) 

the  differential  equation  form  of  the  non-local  dispersion  relation,  equations 
(3.21)  and  (3.22),  reduces  to  the  local  algebraic  dispersion  relation  given  by 
(16)  of  Glassman  and  Sperling  [1983], 
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where 


6 


k^<  T 
K||  B'eo 


(3.33) 


For  spread-F  and  barium-cloud  plasmas,  Glassman  and  Sperling  [1983] 
analyze  the  solution  to  the  local  dispersion  relation,  (3.32),  while  Sperling 
and  Glassman  [1983]  compare  the  solution  to  the  nonlocal  dispersion  relation, 
(3.21),  with  the  solution  to  the  local  dispersion  relation.  A  major  conclusion 
is  that  the  frequency  becomes  increasingly  insensitive  to  the  parallel  wave 
number  as  increases.  Hence,  a  local  flute  estimate  for  the  frequency  of  the 
instability  is  very  reasonable  for  sufficiently  large  k  . 

Although  the  eigenvalue  to  can  be  well  approximated  using  the  flute 
approximation  for  large  k  ,  it  does  not  mean  that  the  eigenfunction,  4',  is 
independent  of  z.  Quite  to  the  contrary,  is  strongly  peaked  in  z  and  is 
confined  to  a  spatial  region  where  the  background  plasma  parameters  do  not  change 
too  substantially  [Sperling  and  Glassman,  1983].  In  fact,  a  local  flute  estimate 
for  the  growth  rate  should  be  expected  to  be  reasonable  when  |k||l_2|  £  1  where  L£ 
is  the  characteristic  length  for  parallel  variations  in  background  plasma 
parameters  [Kadomtsev,  1965]. 


In  the  flute  limit  (i.e.,  k^  =  0),  (3.32)  reduces  to 


[wf(kifz)  -  u0  +  ivjx. 

^(k^z)  -  w0  +  ivc  +  iviAi  “  uyfk^z)  +  iv£ 


which  in  turn  can  be  written  as  a  quadratic 


0  =  [ujf(k^  ,z)]  +  £jW^(k^,z)  +  ^2 


(3.34) 


(3.35) 


In  (3.34)  and  (3.35),  w^(k^,z)  is  the  eigenvalue  as  determined  from  the  flute 
approximation  and  is  a  function  of  kj_  and  z.  In  (3.35) 

C1  =  -u0  +  ivi  +  ivc 


r  .  “O^i 

2  "  TT 


ivTu>i  -  iuovc 


(3.36) 
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The  solution  to  (3.35)  is 


“f  (kj.  *z)  =  *  2 


a  ■  a 


(3.37) 


Comparisons  between  (3.37)  and  numerical  solution  for  w  using  (3.21) 
verify  that  for  a  sufficiently  large  value  of  k 


wf(kl»2m) 


(3.38) 


In  (3.38)  z  is  the  value  of  z  which  satisfies  the  condition  that 
m 


R ( k± , z )  =  Im[ujf(ki,2)]  -  Re[u>f(ki,z)] 


(3.39) 


is  a  maximum,  i .e. , 


3R<kl’Zm) 

3z„ 


aR(ki.zm) 

__ 


(3.40) 


The  ability  of  ojf ( k±  ,zm)  to  closely  approximate  the  actual  eigenvalue  to 
suggests  a  simple  algorithm  for  estimating  an  outer-scale  wave  number  without 
actually  having  to  solve  the  differential  equation,  (3.21).  If  it  is 
hypothesized  that  appreciable  striations  require  Im(u>)  £  |Re(u>) |  ,  then  k^  can  be 
identified  with  the  outer-scale  wave  number.  Equations  (3.39),  (3.40),  and  the 
constraint 


0  ■  R(kt,zm) 


(3.41) 


then  completely  determine  k^.  The  algorithm,  used  to  determine  kt,  is  summarized 
in  Figure  3.1. 

In  Appendix  G  explicit  expressions  are  derived  for  not  only  kt  but  also 
the  perpendicular  wave  number  appropriate  to  zero  growth  rate  (i.e.,  k$). 
Equation  (3.37)  is  used  in  the  evaluation  of  the  various  relevant  expressions. 


o 


Figure  3.1.  Algorithm  for  estimating  the  outer-scale  wave  number 
kt.  The  numbers  designate  equations  in  the  text. 


In  addition  to  being  able  to  make  reasonable  estimates  of  actual 
eigenvalues  in  the  limit  of  large  k^,  it  is  also  possible  to  make  good  estimates 
for  the  spatial  extent  of  the  parallel  density  variations  of  linear  striation 
eigenmodes  (i.e.,  |nej|)  when  k  is  large.  When  the  region  of  growth  is  peaked 
at  high  altitude,  above  the  ionosphere,  then 


I  nel I  88  nfi  5  Kloneo“f(k i  ,z) I 


(3.42) 


with  a  specified  by  (3.17b)  and  K  a  constant.  However,  if  the  region  of  growth  is 
confined  to  lower  altitudes,  within  the  ionosphere,  then 


Vl  1  +  ujn/|af(k,  ,2, 


(3.43) 


The  expressions  for  n^  and  n*2  are  summarized  in  Table  3.2.  If  the  formation  of 
striations  can  be  conceptualized  as  a  mechanism  for  diffusing  late-time  nuclear 
plasmas  then  nfl  and  nf2  provide  estimates  for  the  locations,  along  the 
geomagnetic  field,  where  the  diffusion  occurs. 

3.4  NUMERICAL  RESULTS 

In  order  to  determine  the  eigenvalues  and  eigenmodes  appropriate  to 
late-time,  high-altitude,  nuclear  plasmas,  (3.21)  has  been  explicitly  solved 
using  a  hybrid  col location-Galerk in  method  [Sperling  and  Glassman,  1983]  and 
plasma  profiles  appropriate  to  nuclear  plumes  [r.  E.  Stoeckly  and  R.  W.  Kilb, 
private  communication,  1981;  W.  W.  White,  private  communication,  1982].  The 
parameters  and  geometry  for  the  five  cases  are  summarized  in  Table  3.3.  The 
first  four  cases  differ  only  in  the  values  of  two  parameters:  the  magnitude  of 
the  perpendicular  density-gradient  scale  length,  L^,  and  the  effective  E  x  B 
velocity,  I g^l /vin -  For  Cases  1-4,  ion-gravity  and  ion-curvature  drifts  con¬ 
tribute  to  instability  (i.e.,  gCg  f  0).  Case  5  is  identical  to  Case  1  except 
that  both  the  ion-gravity  and  ion-curvature  drifts  are  artificially  set  equal  to 
zero  (i.e.,  gCg  =  0)  and  so  instability  can  only  be  driven  at  lower  altitudes,  in 
the  ionosphere  where  ion-Pedersen  currents  can  be  appreciable.  Case  5  is 


%*  v  O1  » 


Table  3.3.  Plasma  parameters  for  the  cases  discussed  in  Section  3.4  and 
Appendix  H  .  Subscripts  1  and  2  represent  ions  in  the  unstable  and 
stable  plasmas,  respectively.  For  all  cases  n.  (z)  =  0*5  j  1  + 

tanh [{z  -  ZjJ/Lj  j  +  n-,  nno(z)  =  "no  (Ln/z)7‘  *  and  Li(z)  =  Ui  1 1 

+  exp  [-2  (z  -  ZjJ/Lj]}  jl  +  exp[2  (z  -  ZgJ/Lgl!  .  The  atomic  mass  of 
the  ions  is  assumed  to  be  20.  The  distance  z  is  in  kilometers.  The 
entries  for  YQ  and  Y  assume  that  the  perpendicular  density-gradient 
scale  length  is  L  .  The  altitude,  h,  can  be  related  to  z  by 
h  =  0.92  z. 


cm  ) 
cm-3) 


km/s) 

's2) 


>o  <eV> 


Case  1 

Case  2 

Case  3 

Case  4 

1  x  107 

1  x  107 

1  x  107 

1  x  107 

2.5  x  102 

2.5  x  102 

2.5  x  102 

2.5  x  102 

3.5  x  101 

3.5  x  101 

3.5  x  101 

3.5  x  101 

1  x  104 

1  x  104 

1  x  104 

1  x  104 

4.18  x  108 

4.18  x  108 

4.18  x  108 

4.18  x  108 

3  x  102 

3  x  102 

3  x  102 

3  x  102 

1.5  x  101 

1.5  x  101 

1.0 

1.0 

2  x  103 

2  x  103 

2  x  103 

2  x  103 

1  x  103 

1  x  103 

1  x  103 

1  x  103 

0.3 

0.3 

0.3 

0.3 

0.1 

1.0 

0.1 

1.0 

6.7  x  10-3 

6.7  x  10"3 

6.7  x  10'3 

6.7  x  10"3 

4.3  x  10'3 

4.3  x  10"3 

4.3  x  10'3 

4.3  x  10"3 

1.7  x  10-2 

1.7  x  10-2 

1.7  x  10-2 

1.7  x  10'2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

Case  5 


2.5  x  10‘ 

3.5  x  10J 

1  x  104 
4.18  x  10 

3  x  102 

1.5  x  10] 

2  x  103 
1  x  103 

0.3 

0.1 

6.7  x  10" 


discussed  in  Appendix  H.  The  expressions  and  parameters  shown  in  Table  3.3  are 

based  upon  information  appropriate  to  altitudes  above  100  km.  The  table 

7  -3 

indicates  that  the  plasma  density  is  largest  (i.e.,  ~  1  x  10  cm  )  for 

2 

altitudes  above  ~  2.5  x  10  km  and  that  the  perpendicular  density-gradient  scale 

length  is  smallest  (i.e.,  ~  15  km  for  Cases  1,  2,  and  5  and  —  1  km  for  Cases  3 

2  3 

and  4)  for  altitudes  between  ~2.5  x  10  km  and  2  x  10  km. 

Now  for  intermediate  latitudes  (e.g.,  9  =  46°)  and  an  altitude  of 
3  3 

1  x  10  km,  R  %  7.4  x  10  km  in  (3.1).  If  in  addition  w  ^  ^  5  km/s  and  the 
inverse-square  law  for  gravity  is  used  [ K i 1 b ,  1977],  (3.1)  suggests  that  a 
reasonable  value  for  ) 9Cg  |  is  4.3  x  10”^  km/s^.  This  value  of  1 9Cg I i s  specified 
in  Table  3.3. 

Figures  3.2a,  3.2b,  3.2c,  and  3. 2d  are  graphs  of  w  (s~^)  versus  k.  (km~*) 
for  Cases  1,  2,  3,  and  4,  respectively.  The  perpendicular  wave  number  for 
maximum  growth  rate  varies  between  6  km"^  for  Case  1,  the  most  weakly  driven  of 
the  four  cases,  to  about  30  km“*  for  Case  4,  the  most  strongly  driven  of  the  four 
cases.  The  peak  growth  rate  shifts  to  shorter  perpendicular  wavelengths  as  the 
strength  of  the  driving  currents  increases  relative  to  classical  ambipolar 
diffusion. 

As  kx  is  decreased  from  the  value  for  maximum  growth,  parallel  currents 
act  to  reduce  the  striation  growth  rate.  But  as  k^  becomes  still  smaller,  the 
growth  rate  levels  off  at  a  fraction  of  the  maximum  growth  rate.  This  coincides 
with  the  electromagnetic  effects  becoming  important.  If  only  electrostatic 
terms  are  kept  in  the  solution,  the  growth  rate  continues  to  decline  with 
decreasing  k^,  as  shown  in  Appendix  H. 

As  k^  is  increased  from  the  value  for  maximum  growth,  the  growth  rate  is 
decreased  and  stability  is  eventually  attained  because  of  the  combination  of  ion- 
viscosity  and  classical  ambipolar  diffusion.  The  contribution  of  ion-viscosity 
to  the  growth  rate  is  examined  in  Appendix  H.  As  kj_  increases,  the  real  part 
of  u>  grows  rapidly  and  equals  the  growth  rate  for  kj_  =  k^.  As  discussed  in 
section  3.3,  k^.  is  hypothesized  to  be  the  effective  outer  scale  wave  number 
appropriate  to  the  "freezing"  of  striations.  In  Appendix  G,  "it  is  shown  that 
kt  can  be  appreciably  smaller  than  the  perpendicular  wave  number  appropriate  to 
absolute  stabilization,  k  ,  when  oo*  »  vr.  Thus,  it  is  in  Cases  3  and  4 
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Figure  3.2.  u>  (s'1)  versus  k±  (km-1)  for  (a)  Case  1,  (b)  Case  2, 
(c)  Case  3,  and  (d)  Case  4.  The  solid  and  dashed  lines 
designate  the  growth  rate  and  real  frequency, 
respectively. 


(Figures  3.2c  and  3. 2d),  where  L  is  the  smallest  and  u>*  is  largest  for  specified 
ki,  that  the  effect  of  Re(ui)  has  the  greatest  effect  on  differentiating  between 
kt  and  k$.  The  values  for  kt  generated  using  the  flute  approximation  algorithm 
described  in  section  3.3  are  20,  60,  35,  and  110  km”*  for  Cases  1,  2,  3,  and  4, 
respectively.  This  agrees  well  with  the  more  exact  results  shown  in  Figure  3.2 
(i.e.,  kt  %  20,  50,  35,  and  100  km”*  for  Cases  1,  2,  3,  and  4,  respectively). 

It  was  pointed  out  in  section  3.3  that  as  kj_  becomes  large,  the  solution 
of  the  local  dispersion  relation,  (3.32),  becomes  increasingly  insensitive  to 
k || .  The  density  of  modes  in  the  phase  space  of  effective  parallel  wave  numbers 
is  high  and  so  it  is  not  surprising  that  the  numerical  determination  of  eigen- 
modes  and  eigenvalues  becomes  more  difficult  for  larger  values  of  k^  (especially 
for  k^  &  kt).  For  k^  £  kt,  the  eigenvalues  are  rather  insensitive  to  the 
large-scale  variation  of  plasma  parameters  along  the  geomagnetic  field. 

Figures  3.3a,  3.3b,  3.3c,  and  3.3d  are  graphs  of  lE^I/k^  versus  z  (km) 

for  Cases  1,  2,  3,  and  4,  respectively.  Curves  for  several  different  values  of 

kx  (km”*)  are  plotted  in  each  figure.  As  in  the  nonnuclear  examples  discussed  by 

Sperling  and  Glassman  [1983],  the  striation  eigenmodes  become  more  localized 

along  the  geomagnetic  field  as  the  perpendicular  wave  number  becomes  larger. 

Ion-polarization  currents  are  essential  to  the  localization  at  higher  altitudes 

where  ion-Pedersen  currents  are  negligible  [Sperling,  1983e] .  Evidently  in 

late-time  HANE  plasmas,  localized  disturbances  along  the  geomagnetic  field 

generally  lead  to  a  localized  response  and  so  the  geomagnetic-field  lines  are  not 

equipotentials.  The  peaks  of  the  eigenmodes  shift  to  lower  altitudes  as  the 

perpendicular  wave  number  increases.  This  effect  is  most  evident  in  Figure  3.3a 

(Case  1)  and  in  Figure  3.3c  (Case  3),  the  two  cases  corresponding  to  the  smaller 

values  of  ! gE I / v-jn  in  Table  3.3.  Cases  1  and  3  show  a  tendency  for  collisionless 

instabilities  driven  by  ion-gravity  and  ion-curvature  drifts  at  higher  altitudes 

above  the  neutral  atmosphere.  For  Case  1  (Figure  3.3a)  the  eigenmodes  appear  to 

be  especially  collisionless  and  peaked  at  higher  altitudes  for  k  not  too  near 

1  -1/2 

stabilization.  Now  recall  that  the  characteristic  growth  rates  go  as  L,  and 
for  collisionless  and  collisional  instabilities,  respectively  [cf,  (3.4)  and 
(3.5)].  In  comparing  Figures  3.3a  and  3.3c,  it  is  then  reasonable  that  the  modes 
in  Case  1  (Figure  3.3a)  with  =  15  km  tend  to  be  more  collisionless  and  extend 
to  higher  altitudes  than  in  Case  3  (Figure  3.3c)  with  1.^  =  1  km.  Similarly,  the 


lEiJ/k,  versus  z  (km)  for  various  k.  (km"  )  and  for 
(a]yCase  1,  (b)  Case  2,  (c)  Case  3,  and  (d)  Case  4. 

Each  of  the  curves  are  in  normalized  dimensionless 


eigenmodes  shown  in  Figure  3.3b  (Case  2)  extend  to  higher  altitudes  than  in 
Figure  3.3d  (Case  4).  The  tendency  for  the  peak  of  the  eigenmodes  to  shift  to 
lower  altitudes  as  perpendicular  wave  number  increases  means  that  striations  at 
higher  altitudes  will  tend  to  have  larger  perpendicular  scale  sizes  than 
striations  at  lower  altitudes.  It  also  means  that  for  the  case  when  both  and 
gCg  are  both  finite,  the  eigenmodes  become  more  collisional  and  tend  to  peak  at 
lower  altitudes  as  kj_  increases  and  approaches  kt.  Hence,  in  Appendix  G, 
formulas  for  kt  are  reasonably  derived  assuming  the  collisional  limit  when  ion- 
neutral  collisions  are  important.  The  high  altitude  peaking  of  Case  1  and  its 
relation  to  ion-Pedersen  versus  ion-gravity  and  ion-curvature  drifts  is  explored 
further  in  Appendix  H. 

The  solutions  to  (3.35)  [i . e - ,  ( k^ , z )  (s-1)]  are  plotted  as  a  function 

of  z  (km)  in  Figures  3.4a,  3.4b>  3.4c,  and  3.4d,  for  Cases  1,  2,  3,  and  4,  and  k 

_1  1 
values  of  10,  30,  30,  and  70  km  ,  respectively.  The  k±  values  are  chosen  to  be 

larger  than  that  appropriate  to  maximum  growth  but  smaller  than  k^.  The  values 

for  zm  can  be  calculated  from  (3.39)  and  (3.40)  and  are  approximately  690,  280, 

280,  and  260  km  for  Cases  1,  2,  3,  and  4,  respectively.  The  values  of  zm  roughly 

match  the  appropriate  values  for  z  where  lE^j/^  are  maximum  in  Figures  3.3a-d. 

The  approximate  eigenvalues  [i .  e . ,  u>f  (ki,zm)  (s”^)j,  deduced  from 
(3.35),  (3.39),  and  (3.40)  are  plotted  as  dashed  lines  versus  (km-1)  in 
Figures  3.5a,  3.5b,  3.5c,  and  3.5d,  for  Cases  1,  2,  3,  and  4,  respectively.  The 
solid  lines  in  Figures  3.5a-d  are  gj(s-1}  versus  k^  (km-1),  replotted  from 
Figures  3.2a-d.  The  curves  for  uy  (k^,zm)  and  u>  become  especially  close  for 
k  &  kt.  This  demonstrates  that  the  flute  approximation  yields  accurate 
estimates  for  the  eigenvalues  vhen  k^  becomes  large  enough.  For  small  values  of 
k^ ,  eigenmodes  sample  large  variations  in  background  plasma  parameters. 
Im  [co^  ( k^  ,zm)]  then  overestimates  the  growth  rate  as  parallel  electron  currents 
act  to  short-circuit  the  perpendicular  electric  fields  in  eigenmodes. 

Figures  3.6a,  3.6b,  3.6c,  and  3.6d  are  plots  of  |ngJ  ,  n^,  and  n^ 

versus  z  (km)  for  Cases  1,  2,  3,  and  4,  respectively.  Expressions  for  n^  and 

n^2  are  given  by  (3.42)  and  (3.43),  respectively.  The  empirical  estimates,  n^ 

for  Case  1  and  n^  for  Cases  2,  3,  and  4,  give  reasonable  approximations  to  |neJ 

particularly  at  low  altitudes  where  |nej|,  n^,  and  n ^  mirror  the  spatial 

variation  in  n  . 
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Figure  3.4.  cof  (k± , zm)  versus  z  (km)  for  (a)  Case  1,  (b)  Case  2, 
(c)  Case  3,  and  (d)  Case  4.  Specific  values  of  kj_ 
(km-1)  are  assumed.  The  solid  and  dashed  lines 
designate  the  growth  rate  and  real  frequency, 
respectively. 
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Figure  3.6. 


I nel  *  nfl>  and  nf2  versus  z  (km)  for  (a)  Case  1, 
(b)  Case  2,  (c)  Case  3,  and  (d)  Case  4.  The  solid  line 
is  | n  J  as  determined  from  (3.17a)  and  the  solution  of 
(3.21J.  The  dashed  line  is  (a)  nfl  and  (b)-(d)  nf2. 
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3.5  SUMMARY  AND  CONCLUDING  REMARKS 

A  second-order  linear  differential  equation  has  been  derived  for 
determining  the  linear  eigenmodes  of  striation  instabilities  along  the 
geomagnetic  field  and  has  been  solved  for  parameters  and  geometry  appropriate  to 
late-time  HANE  plasmas,  as  described  by  MELT.  Eigenvalues  (i.e.,  frequencies) 
are  determined  along  with  the  eigenmodes.  The  differential  equation  includes  the 
contribution  of  two  physical  effects  which  are  essential  for  correctly 
determining  the  eigenvalues  and  the  localization  of  the  eigenmodes.  These 
effects  are  ion-polarization  currents  and  inductive  electric  fields.  Ion- 
polarization  currents  are  especially  important  at  higher  altitudes  where  ion- 
Pedersen  currents  are  negligible.  Ion-polarization  currents  are  a  mechanism  for 
short-circuiting  voltage  differences  between  adjacent  geomagnetic-field  lines. 
Inductive  effects  are  required  in  order  to  properly  take  into  account  the  fact 
that  information  is  propagated  along  the  geomagnetic  field  at  a  finite  rate  which 
depends  upon  the  Alfven  speed. 

It  is  demonstrated  that  striations  tend  to  develop  with  larger  perpen¬ 
dicular  scale  sizes  at  higher  altitudes  than  at  lower  altitudes.  The  results 
also  demonstrate  that,  in  late-time  HANE  plasmas,  localized  disturbances  along 
the  geomagnetic  field  generally  lead  to  a  localized  response  and  so  the 
geomagnetic-field  lines  are  not  necessarily  equipotentials.  Electromagnetic 
effects  are  shown  to  be  important  for  modes  with  smaller  perpendicular  wave 
numbers  (i.e.,  Si  km”*)  and  ion-viscosity  is  important  for  modes  with  larger 
perpendicular  wave  numbers  (i  •  6  •  *  &  4  km-*).  Analytic  criteria  are  given  for  the 
perpendicular  wave  numbers  of  the  instability  at  which  the  real  part  of  the 
frequency  is  equal  to  the  growth  rate  and  at  which  there  is  complete  stabiliza¬ 
tion.  It  is  argued  that  the  equating  of  real  part  of  the  frequency  and  the  growth 
rate  gives  a  meaningful  estimate  of  the  perpendicular  wave  number  at  which 
striation  "freezing"  occurs.  An  appropriate  algorithm,  suitable  for  implemen¬ 
tation  in  modeling  codes,  is  given  for  determining  the  perpendicular  wave  number 
at  "freezing."  Simple  empirical  estimates  for  the  parallel  spatial  extent  of  the 
striation  density  fluctuations  are  also  described. 


REFERENCES 


Alfven,  H.,  Cosmic  Plasma,  D.  Reidel,  Boston,  1981. 

Allis,  W.  P.,  S.  J.  Buchsbaum,  and  A.  Bers,  Waves  in  Anisotropic  Plasmas,  MIT 
Press,  Cambridge,  Massachusetts,  1963. 

Block,  L.  P.,  Potential  double  layers  in  the  ionosphere.  Cosmic  Electrodynam. ,  3, 
349,  1972. 

Boyd,  T.  J.  M.,  and  J.  S.  Sanderson,  Plasma  Dynamics,  Barnes  and  Noble,  New  York, 
1969,  pp.  51-53. 

Braginskii,  S.  I.,  Transport  processes  in  plasmas,  in  Reviews  of  Plasma  Physics, 
Vol.  1,  edited  by  M.  A.  Leontovich,  Consultants  Bureau,  New  York,  1965. 

Glassman,  A.  J.,  and  J.  L.  Sperling,  Electromagnetic  theory  of  collisional 
interchange  instabilities,  J.  Geophys.  Res.,  88,  10091,  1983. 

Goldman,  S.  R.,  L.  Baker,  S.  L.  Ossakow,  and  A.  J.  Scannapieco,  Striation 
formation  associated  with  barium  clouds  in  an  inhomogeneous  ionosphere, 
J.  Geophys.  Res.,  81,  5097,  1976. 

Goldman,  S.  R.,  and  J.  L.  Sperling,  Electric  field  limitations  of  drift  dissi¬ 
pative  wave  generation  in  ionospheric  structures,  J.  Geophys.  Res.,  87,  254, 
1982. 

Hudson,  M.  K.,  and  C.  F.  Kennel,  Linear  theory  of  equatorial  spread  F, 
J.  Geopnys.  Res.,  80,  4581,  1975. 

Jackson,  J.  0.,  Classical  Electrodynamics,  John  Wiley  and  Sons,  New  York,  1962, 
ch.  7. 

Kadomtsev,  B.  B.,  Plasma  Turbulence,  Academic  Press,  New  York,  1965. 

Kilb,  R.  W.,  Striation  formation,  in  Physics  of  High-Altitude  Nuclear  Bursts, 
Defense  Nuclear  Agency  Rep.  0NA  4501F,  Defense  Nuclear  Agency,  Washington, 
D.C.,  1977. 

Knapp,  W.  S.,  and  K.  Schwartz,  Aids  for  the  study  of  electromagnetic  blackout. 
Defense  Nuclear  Agency  Rep.  DNA  3499H,  Defense  Nuclear  Agency,  Washington, 
D.C.,  1975,  p.  8-5. 

Landau,  L.  D.,  and  E.  M.  Lifshitz,  Electrodynamics  of  Continuous  Media,  Pergamon, 
New  York,  1960. 

Linson,  L.  M.,  and  J.  B.  Workman,  Formation  of  striations  in  ionospheric  clouds, 
J.  Geophys.  Res.,  75,  3211,  1970. 


Linson,  L.  M.,  and  G.  Meltz,  Theory  of  ion  cloud  dynamics  and  morphology,  in 
Analysis  of  Barium  Clouds,  Rep.  RADC-TR-72-736,  Ch.  5,  Vol.  1,  Avco  Everett 
Research  Laboratory,  Everett,  Massachusetts,  1972. 

Linson,  L.  M.,  and  D.  C.  Baxter,  Ion  cloud  modelling.  Defense  Nuclear  Agency  Rep. 
DNA  4455F,  Defense  Nuclear  Agency,  Washington,  D.C.,  1977. 

Longmire,  C.  L.,  Elementary  Plasma  Physics,  Interscience,  New  York,  1963. 

McDonald,  B.  E.,  M.  J.  Keskinen,  S.  L.  Ossakow,  and  S.  T.  Zalesak,  Computer 
simulation  of  gradient  drift  instability  processes  in  Operation  Avefria, 
J.  Geophys.  Res.,  85,  2143,  1980. 

McDonald,  B.  E.,  S.  L.  Ossakow,  S.  T.  Zalesak,  and  N.  J.  Zabusky,  Scale  sizes  and 
lifetimes  of  F  region  plasma  cloud  striations  as  determined  by  the  condition 
of  marginal  stability,  J.  Geophys.  Res.,  86,  5775,  1981. 

Mende,  S.  B.,  Morphology  of  •  the  magnetospheric  barium  release,  J.  Geophys.  Res., 
78,  5751,  1973. 

Ossakow,  3.  L.,  S.  T.  Zalesak,  B.  E.  McDonald,  and  P.  K.  Chaturvedi,  Nonlinear 
equatorial  spread  F:  Dependence  on  altitude  of  the  F  peak  and  bottomside 
background  density  gradient  scale  length,  J.  Geophys.  Res.,  84,  17,  1979. 

Ott,  E.,  Theory  of .  Rayleigh-Taylor  bubbles  in  the  equatorial  ionosphere, 
J.  Geophys.  Res.,  83,  2066,  1978. 

Overman,  E.  A.,  II,  N.  J.  Zabusky,  and  S.  L.  Ossakow,  Ionospheric  plasma  cloud 
dynamics  via  regularized  contour  dynamics,  I.  Stability  and  nonlinear 
evolution  of  one-contour  models,  Phys.  Fluids,  26,  1139,  1983. 

Paoloni,  F.  J.,  Boundary  effects  on  M  =  0,  ±1  Alfven  waves  in  a  cylindrical, 
collisionless  plasma,  Phys.  Fluids,  18,  140,  1975. 

Papadopoulos,  K.,  A  review  of  anomalous  resistivity  for  the  ionosphere.  Rev. 
Geophys.  Space  Phys.,  15,  113,  1977. 

Perkins,  F.  W.,  N.  J.  Zabusky,  and  J.  H.  Doles,  III,  Deformation  and  striation  of 
plasma  clouds  in  the  ionosphere,  J.  Geophys.  Res.,  78,  697,  1973. 

Prettie,  C.  W.,  Results  of  Places  data  analysis.  Defense  Nuclear  Agency  Rep.  DNA- 
TR-81-270,  Defense  Nuclear  Agency,  Washington,  D.C.,  1982. 

Prettie,  C.  W.,  S.  Y.  F.  Chu,  and  J.  B.  Workman,  The  microstructure  theory  for 
three-dimensional  plasma  transport.  Defense  Nuclear  Agency  Rep.  DNA-TR-81- 
270,  Defense  Nuclear  Agency,  Washington,  D.C.,  1983. 

Scannapieco,  A.  J.,  and  S.  L.  Ossakow,  Nonlinear  equatorial  spread  F,  Geophys . 
Res.  Lett..  3,  451,  1976. 


Scholer,  M.,  On  the  motion  of  artificial  ion  clouds  in  the  magnetosphere.  Planet. 
Space  Sci.,  18,  977,  1970. 

Simon,  A.,  Instability  of  a  partially  ionized  plasma  in  crossed  electric  and 
magnetic  fields,  Phys.  Fluids,  6,  382,*  1963. 

Sperling,  J.  L.,  Inhomogeneity  effects  on  wave  properties  in  Elmo  bumpy  torus 
scale,  Phys.  Rev.  A,  26,  1618,  1982a. 

Sperling,  J.  L.,  Role  of  Rayleigh-Taylor  instabilities  on  prompt  striation 
evolution,  J.  Geophys.  Res.,  87,  10514,  1982b. 

Sperling,  J.  L.,  and  A.  J.  Glassman,  Striation  eigenmodes  along  the  geomagnetic 
field  and  eigenvalues  in  the  limit  of  strong  ion-neutral  collisions,  JAYCOR 
Rep.  J530-83-135,  JAYCOR,  San  Diego,  California,  1983. 

Sperling,  J.  L.,  The  role  of  ion-neutral  diffusion  on  collisional  structuring 
instabilities,  JAYCOR  Rep.  J530-83-144,  JAYCOR,  San  Diego,  California,  1983a. 

Sperling,  J.  L.,  Ion-Pedersen  drift  and  parallel  electric  field  effects  on  plasma 
jetting,  J.  Geophys.  Res.,  88,  7095,  1983b. 

Sperling,  J.  L.,  On  the  effect  of  finite,  field-aligned  plasma  length  on  a  loss- 
cone  instability,  J.  Geophys.  Res.,  88,  927,  1983c. 

Sperling,  J.  L.,  Finite  parallel  wavelengths  and  ionospheric  structuring, 
J.  Geophys.  Res.,  88,  4075,  1983d. 

Sperling,  J.  L.,  The  parallel  evanescence  of  striation  parameters,  JAYCOR  Rep. 
J530-83-119,  JAYCOR,  San  Diego,  California,  1983e. 

Spiegel,  M.  R.,  Schaum's  Outline  of  Theory  and  Problems  of  Laplace  Transforms, 
McGraw-Hill,  New  York,  1965. 

Stix,  T.  H.,  The  Theory  of  Plasma  Waves,  McGraw-Hill,  New  York,  1962. 

Stratton,  J.  A.,  Electromagnetic  Theory,  McGraw-Hill,  New  York,  1941. 

Vesecky,  J.,  J.  Chamberlain,  J.  Cornwall,  D.  Hammer,  and  F.  W.  Perkins, 
Irregularities  in  ionospheric  plasma  cloud:  Their  evaluation  and  effect  on 
radio  communications,  JASON  Tech.  Rep.,  JSR-80-15,  SRI  International, 
Arlington,  Virginia,  1980. 

Volk,  H.  J.,  and  G.  Haerendel,  Striations  in  ionospheric  ion  clouds,  1, 
J.  Geophys.  Res.,  76,  4541,  1971. 

Zalesak,  S.  T.,  and  S.  L.  Ossakow,  Nonlinear  equatorial  spread  F:  Spatially 
large  bubbles  resulting  from  horizontal  scale  initial  perturbations, 
J.  Geophys.  Res.,  85,  2131,  1980. 


APPENDIX  A 

MATHEMATICAL  DESCRIPTION  OF  CONDUCTIVITIES 


This  appendix  specifies  the  mathematical  forms  for  oy,  o  ,  and  o  in 
_  _  _  x  y  z 

(1.3)  and  a  ,  o  ,  and  o  in  (1.6a-b). 
x  y  z 

The  perpendicular  currents,  driven  by  the  perpendicular  electrostatic 
fields,  are  considered  to  be  a  consequence  of  ion-neutral  collisions  and  temporal 
variations  in  the  electrostatic  fields.  Hence,  the  perpendicular  currents  are  of 
the  Pedersen  or  polarization  types.  The  parallel  current  is  assumed  to  be 
controlled  by  the  parallel  electrostatic  field  and  electron  resistivity.  Hence 
[Sperling,  1983b], 


In  (A2)  the  following  symbols  are  used:  e  (proton  charge),  n  (density), 
c  (speed  of  light),  B  (magnetic  field  strength),  v^n  (ion-neutral  collision 
frequency),  ( ion-gyrofrequency) ,  mg  (electron  mass),  and  vg  (effective 
collision  frequency  for  parallel  electron  motion).  In  writing  (A2)  it  has  been 
assumed  that  vip,  |3/3t|  «  and  that  ions  have  the  proton  charge.  Relevant 
mathematical  forms  for  v..^  and  vg  are  [Braginskii,  1965;  Ki lb,  1977] 
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23.4  -  1.15  log1Q(n)  +  3.45  1og10(Te) 
3.5  x  105 


In  (A3)  the  following  symbols  are  used:  m^  (ion  mass),  nn  (neutral  density), 
and  Te  (electron  temperature).  All  parameters  are  in  cgs  units  except  for  the 
electron  temperature  which  is  in  units  of  eV. 


It  is  evident  from  (A2)  that  Oj_  is  affected  by  temporal  variations  in 

the  electrostatic  fields  but  oz  is  unaffected.  The  application  of  the  Laplace 

transform  (1.5a)  to  (1.4a-b)  implies  the  following  functional  forms  for  o  ,  o  , 

_  *  y 

and  oz  in  (1.6a-b) : 
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APPENDIX  B 

TOTAL  ELECTROSTATIC  POTENTIAL  FOR  AN  APPLIED  ELECTROSTATIC  FORCE 


This  appendix  uses  (1.21a-b)  to  derive  the  total  electrostatic  potential 
for  the  case  when  the  external  applied  force  is  completely  electrostatic. 
Specifically, 
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The  electrostatic  potential  for  the  applied  electrostatic  field  is 
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ox 
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( B2) 


Using  the  Laplace  transformation,  (1.5a),  permits  the  total  electrostatic 
potential  to  be  written  as 
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(B3) 


The  combination  of  (1.21a),  (1.21b),  (Bl),  and  (B2)  gives 
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( B4b) 


Equation  (B4b)  agrees  with  (35)  and  (38)  on  pages  212-213  of  Stratton  [l94l] 
and  the  result  of  problem  6  on  page  62  of  Landau  and  Lifshitz  [i960]. 

Equations  (B2)  and  (B4a)  indicate  that  the  electric  field  strength 
in  the  interior  of  the  ellipsoid  is  reduced  (increased)  relative  to  the  value  far 
from  the  ellipsoid  if  the  appropriate  conductivity  inside  the  ellipsoid  is  larger 
(smaller)  than  the  appropriate  conductivity  outside  the  ellipsoid.  For  example, 
V  >(<)  °i>  implies  that  the  x-directed  electric  field  is  smaller  (larger) 
inside  the  ellipsoid  than  far  from  the  ellipsoid. 
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The  final-value  theorem  for  Laplace  transformations  indicates  that  the  steady 
state  value  for  the  parallel  electric  field  on  the  surface  of  the  ellipsoid  is 
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On  the  surface  of  the  ellipsoid 


<2  +  y2  *  a2(l  -  z2/c2) 


2  2 

Hence,  for  y  =  0  and  o^v^Cq/o  ag  <K  1> 
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The  magnitude  of  the  right-hand  side  has  a  maximum  value  for 
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To  within  a  constant  (F6)  is  in  agreement  with  (1.21a)  and  the  appropriate 
asymptotic  form  of  entries  3  and  4  in  Table  1.2. 

0  5 

If  b^Jo  )  "  >:>  c0*  the  field-line  integrated  Pedersen  conductivities 
of  the  plasma  cloud  and  the  ambient  plasma  are  determined  from  (F6)  to  be 
co(°L<  )°*5  and  bo(0z>)°*5ai>’  respectively.  Hence,  the  concept  of  field-line 
integrated  Pedersen  conductivities  appears  to  be  approximately  applicable,  to 
within  factors  of  order  unity,  for  at,  ct  «  bt  or  ct  «  at  =  bt. 
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APPENDIX  F 

INDUCED  ELECTRIC  FIELDS  AND  FIELD-LINE  INTEGRATION 


It  is  evident  from  (1.18),  (1.22a),  (1.22b),  and  entries  (3)-(5)  of 
Table  1.2,  that  the  depolarization  factors  are  complicated  functions  of  geometry 
when  the  semiprincipal  axes  (i.e.,  a^,  bt,  and  ct)  are  of  comparable  magnitude. 
Under  such  circumstances,  field-line  integration  would  appear  to  be 
inappropriate  for  determining  induced  electric  fields  inside  plasma  clouds. 

Consider  the  situation  of  j  >  j-  0,  j  >  =  j  >  =  0,  and  c. ,  at  «  b.  or 
ct  «  a^  =  bt.  Overall  quasineutrality  requires  that  net  current  flow  into  and 
out  of  the  ellipsoid  be  nondivergent.  The  total  x-directed  current  inside  the 
plasma  cloud  is 

ll  *  ai<Ex<1lb0c0  +  J‘ox<  •  (pl) 


Equation  (FI)  uses  the  fact  that  the  total  surface  area  presented  to  the 
x-directed  current-density  is  ~  nbgCg.  The  total  x-directed  current  flowing 
from  the  external  plasma  directly  into  the  plasma  cloud  is 

'z  *  °x>Ex>"b0c0  +  Jox>  •  (F2) 

To  maintain  closed  current  loops,  there  is  a  return  current  in  the  external 

plasma  which  is  present  above  and  below  the  plasma  cloud  with  respect  to  position 

along  the  geomagnetic  field  lines.  Considerations  similar  to  those  leading  to 

1/2 

( 1 . 29 ) -( 1 .31 )  indicate  that  the  return  current  extends  a  distance  ~~  bt/az>^ 
along  the  magnetic  field  line.  Hence,  the  return  current  is 
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Total  conservation  of  current  requires 


less  likely  to  occur  for  plasma  clouds  with  transverse  scale  size  approaching  L  . 

When  considering  the  formation  of  striations  L  can  be  equated  with  an  inverse 

-1  n  4 

perpendicular  wave  number,  k  .  If  cE  /B  *  1  x  10  cm/s,  B  =  0.3  G,  and 

i  y  3 

T.  =  0.1  eV,  then  Ln  %  k  =  3.3  x  10  cm  and  the  corresponding  transverse 
1  n  y  4 

wavelength  is  =  2n/ky  »  2.1  x  10  cm.  This  value  for  is  comparable  to  the 
observed  scale  size  for  "frozen"  striations  [Linson  and  Meltz,  1972]. 


Equation  { E5 )  assumes  that  electric  fields  are  not  induced.  The  density 
profile  implied  by  (E4)  is 
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with  C  a  constant  of  integration.  If  v  is  an  ion-Pedersen  velocity  then 

,  U  A 


V,  ffox 

'ox  ‘  Q.  B 


and  (E6)  reduces  to 
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Equations  (E8)  and  (E9)  suggest  that  Ln  is  an  important  scale  size  in 
the  dynamics  of  plasma  clouds  when  ion-neutral  collisions  are  important.  If  a 
plasma  cloud  has  a  much  larger  dimension  than  Ln  across  the  magnetic  field,  then 
electric  fields  rather  than  ion-neutral  diffusion  can  be  expected  to  be  the 
primary  mechanism  contributing  to  the  maintenance  of  overall  nondivergent 
current  densities.  However,  as  the  perpendicular  scale  size  approaches  Ln,  then 
ion-neutral  diffusion  can  more  likely  be  expected  to  contribute  to  the 
maintenance  of  overall  nondivergent  current  densities. 

Striations  are  observed  to  elongate  outward  from  the  backside  of  barium 
clouds.  Electric  fields  are  required  for  this  bulk  plasma  motion,  but  it  appears 
that  when  the  perpendicular  scale  of  a  plasma  cloud  approaches  Lp  the  tendency  to 
induce  electric  fields  can  be  lessened.  Hence,  the  formation  of  striations  is 
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APPENDIX  E 

DIVERGENT  CURRENT  DENSITIES  AND  ION-NEUTRAL  DIFFUSION 


Figure  D1  shows  a  one-dimensional  plasma  cloud  bounded  at  x  =  ±xQ. 
There. is  a  density  gradient  inside  the  plasma  cloud.  It  is  assumed  that  the 
plasma  external  to  the  cloud  is  uniform  and  symmetric  about  x  =  0.  An  applied 
external  force  results  in  current  densities,  j.v>. 

The  ion-diffusive  current  for  the  geometry  of  Figure  D1  is 


j 


ox 


dn 

-e  Din<  ar 


(El) 


with  the  coefficient  of  ion-neutral  diffusion  given  by 


v.  p? 
in<  i< 


J 


(E2) 


with 


kBTi< 

mi<ni< 


(E3) 


and  kg  the  Boltzmann  constant. 

Consider  a  drift  velocity,  v  ,  which  is  assumed  to  be  independent  of  x. 

U  A 

If 


^ox  "  envox  ’ 

and  if  j  f  j  then  overall  nondivergent  current  densities  away  from 

UA  ^  UA ^ 

| x |  =  xQ  require 


0 


dn. 


<n<  -  n>>vox  -  Din<  dF 


(E5) 
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(Viscosity  is  not  important  if  the  time  rate  of  change  of  Ex  is  larger  than  the 
ion-ion  collision  frequency.)  When  (Oil)  is  satisfied,  the  polarization 
current  density,  rather  than  the  viscous  current  density,  is  the  main  contributor 
to  overall  current  flow  within  the  bulk  of  the  plasma  cloud.  If 
n<  =  1  x  10^  cm"^,  =  0.1  eV,  =  20  n^,  BQ  =  0.03  Q,  and  xQ  =  2  x  10^  cm, 
then  )<<Xo|  «  1  is  satisfied  for  |s|  »  7.2  x  10”^  s"^.  Provided  <$^  >  1,  the 
viscous  force  increases  with  altitude  if  all  other  parameters  are  kept  the  same 
and  so  the  parameter  s  must  be  correspondingly  larger  for  (Dll)  to  be 
satisfied.  At  sufficiently  high  altitudes  S<  <  1  and  the  viscous  force  does  not 
depend  upon  the  magnetic  field  strength  within  the  plasma  cloud. 
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As  a  quantitative  example,  consider  a  barium  cloud,  at  an  altitude  of 
200  km,  with  n<  =  lx  10^  cm"^,  n  =  1.2  x  10^  cm"^,  T.>.  =  0.1  eV, 

Q  T  5  '  < 

np  *  8.3  x  10  cm”  ,  and  m^  =  20  mp,  B  =  0.4  G,  and  xQ  =  2  x  10  cm.  At  steady 

state,  (08)  and  (09)  indicate  that  <<  ~  4.5  x  10"^  cm"^  and 

2  1  3 

^  =  5.5  x  10  cm  and  so  the  boundary  layer  extends  ~  2.2  x  10  cm  into  the 

plasma  cloud.  The  values  of  <<  and  <>  are  sufficiently  large  (i.e.,  1 

and  <>Pi>  »  1  with  Pi  the  ion-gyroradius)  so  that  finite  gyroradius  effects,  not 

included  in  the  present  analysis,  would  likely  impact  actual  boundary  layers.  In 

any  case,  it  appears  that  viscosity  permits  the  existence  of  plasma  structures 

with  large  perpendicular  gradients  at  an  altitude  of  200  km.  However,  since  Ex 

varies  near  the  edge  of  the  plasma  it  follows  that  the  E  x  B  velocity  must  also 

vary.  Hence,  viscosity  is  an  additional  mechanism,  besides  diffusion  and 

Pedersen  convection  [Linson  and  Meltz,  1972;  Prettie,  1982],  for  ablating  plasma 

clouds  and  thereby  transforming  plasma  rods  into  plasma  sheets  [Linson  and  Meltz, 

1972]. 

As  altitude  increases,  (D9)  indicates  that  <>  0  since  v.  -*•  0,  and 

z  in 

so  the  boundary  layers  become  thick  at  sufficiently  high  altitudes  and  steady 
state.  If  m^  =  m^<t  =  T^,  6^  »  1,  and  n<  »  n>,  then  (02)  and  (D10) 
indicate  that  <<  «  and  the  viscous  current  density  specified  by  (D3)  is  a 
main  contributor  to  current  flow  throughout  the  plasma  cloud,  not  only  at  the 
edge  of  the  cloud.  Since  ablation  would  appear  to  act  over  a  large  distance 
inside  a  plasma  cloud  if  <<  xQ  <  1,  it  appears  that  large  perpendicular  gradients 
are  then  not  permitted. 

Viscosity  allows  large  perpendicular  gradients  at  the  boundaries  of 
plasma  clouds  if 


I  <<xo  I  .  I  <>xo  I  >:>  1  •  (°10) 

If  the  parameter  s  in  (D6)  is  considered  to  be  representative  of  the  time  rate 
of  change  of  Ex,  then  even  if  v^n  0  sharp  perpendicular  gradients  are  possible 
if 


svii<mi<  _ 

SVii>mi  , 

d<kBTi<  0 

» 

Wi>  0 

»  1 


(Dll) 
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Schematic  illustration  of  viscous  boundary  layers  near  x  *  xQ  for 

|<<|  «  !<>!  ( - )  and  |  <>|  »  |<<|  ( - ).  The  location  of  the 

plasma  cloud  is  |x|  <  xQ  as  shown  in  Figure  Dl. 


For  both  -►  0  and  <»,  the  viscous  force  goes  to  zero. 
(D5a)  and  (D5b)  then  reduce  to 


for  | x |  >  x( 
for  I x |  <  x. 


Equations 


(D7a) 


(D7b) 


Evidently,  the  transition  of  across  x  =  ±xQ  appropriately  occurs  over  an 
infinitesimal  distance. 

For  finite  values  of  v^,  ( 05a )  and  (D5b)  indicate  the  existence  of 
viscous  boundary  layers  with  finite  thickness  at  x  =  ±xQ.  These  boundary  layers 
are  illustrated  in  Figure  D2.  Two  different  situations  arise  depending  on  the 
relative  magnitudes  of  ic,  and  <>.  When  !*<!  *  I  Si  the  electric  field 
appropriate  to  |x|  »  xQ  extends  a  distance  of  approximately  [Re^)]”*  into  the 
plasma  cloud.  Conversely,  when  |k,|  «  |<J  the  electric  field  appropriate  to 
lx|  «  xQ  extends  a  distance  of  approximately  [ReC^)]"  outside  the  plasma 
cloud. 

In  general,  the  temporal  variation  of  Ex  can  be  determined  by  using  the 

inverse  Laplace  transformation,  (1.5b).  Steady-state  values  of  Ex  are 

determined  from  lim  (sEv).  If  j.v  >  are  both  independent  of  time  then 
S  -►  U 


K,  < 

sinh(K<xQ)exp[>c>(x0  -  I x| )]  for  |x|  >  xQ  (D8a) 


lim  E  = 
t  -►  ® 


K1  1 


cosh(<<x) 


for  | x I  <  xQ  (D8b) 


.  |  aK(W  Jox<)  K  f<  sinh(K  x  )  .  CBh(K  X  ) 
1  c  en<  vin<  2  S  •=  0  ■=  0 


*  '  1“  kB  T1  L 


Using  the  Laplace  transformation,  (1.5a),  and  the  perpendicular 
conductivity  specified  by  (A4),  the  requirement,  that  total  current  be  free  of 
divergence,  can  be  written  as 


K  «  j. 


c£  /n  d_\  _  enc  , 
g2  dx  \  1  dx/  Bfli  1  i 


+  s)  E 
in  '  x 


In  (D4),  K  does  not  depend  upon  x  but  may  depend  upon  the  Laplace  transform 
variable,  s. 

A  physically  reasonable  solution  to  (D4)  is  obtained  by  requiring  both 
Ex  and  dEx/dx  to  be  continuous  across  x  *  ±x0.  This  implies  that  the  viscous 
forces  and  the  related  current  are  finite  at  the  boundaries  of  the  cloud. 
Furthermore,  if  E  -►  0  for  |x|  -»-®,  then  it  follows  that  K  =  Jox>* 

The  appropriate  solution  to  (D4)  is 


—  ~  sinhftc,.  x0)exp[<>(xQ  -  |x|)]  for  |x|  >  xQ  (D5a) 


k2  S 


K1 


cosh(<<  x' 


for  1 x 1  <  xQ  (D5b) 


*1  =  f  ^elXn7+°^r  ’  *2  =  ^  sinh(Vo)  +  cosh(<<Xo)  , 


<vin  +  s)viimi  ' 
d  kB  Ti  * 


The  real  part  of  <>  is  chosen  to  be  greater  than  zero. 


Figure  01.  The  one-dimensional  plasma  cloud  used  to  calculate  the  effect  of 
viscosity  and  ion-neutral  diffusion  on  the  creation  of  induced 
electric  fields  resulting  from  currents  traversing  through  an 
inhomogeneous  plasma.  For  the  calculation  of  ion-viscous  effects  the 

density  profile  is  flat  inside  the  plasma  cloud  ( - ),  but  for  the 

calculation  of  ion-neutral-diffusive  effects  a  density  gradient  is 
permitted  inside  the  plasma  cloud  ( - ). 
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APPENDIX  D 

INDUCED  ELECTROSTATIC  FIELDS  IN  A  VISCOUS  MEDIUM 


To  illustrate  the  impact  of  ion  v  icosity  on  the  formation  of  induced 

electrostatic  fields,  the  simplified  one-dimensional  geometry  of  Figure  D1  is 

used.  The  geometry  is  symmetric  about  x  =  0.  A  uniform  plasma  cloud  is  present 

for  0  5  | x |  <  Xq.  The  plasma  external  to  the  cloud  is  also  uniform.  An  applied 

external  force  results  in  current  densities,  jQx^. 

If  ^ox>  ^  Jox<  an  x"directecl  electric  field,  Ex,  is  created.  This 

electric  field  results  in  a  y-directed  plasma  velocity,  v  *  -cEy/B,  and  a 

y  * 

y-directed  viscous  force  per  unit  volume,  i.e.. 


with  [Braginskii,  1965] 


n 


1  = 


nkBTid 


n 


1.26  +  2.23 
52  +  4.036  +  2.33 


> 


vii  = 


23.4  -  1.151og1Q(n)  +  3.451o9l(?(Te) 


3  x  10' 


/2mo\1/2  n 

(; if) 


(D2) 


T^  (ion  temperature),  (ion-ion 


In  (D2),  the  following  symbols  are  used: 

collision  frequency),  and  mp  (proton  mass).  All  quantities  in  (1.3)  are  in  cgs 

units  except  for  temperature  which  is  in  units  of  eV.  The  Boltzmann  constant  k 

-12 

has  the  value  1.6  x  10  ergs/eV.  The  force  specified  by  (Dl)  implies  an 
x-directed  current  density, 


B 


j 


vx 


(D3) 
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APPENDIX  G 

ASYMPTOTIC  FORMS  OF  THE  FLUTE  SOLUTION  AND  APPLICATIONS 


The  quantitative  solution  of  (3.21)  indicates  that  the  flute  approxima¬ 
tion  can  provide  accurate  estimates  for  the  eigenvalues  (i.e.,  frequencies)  for 
large  perpendicular  wave  numbers.  In  this  appendix  (3.36)  and  (3.37)  are  used  to 
derive  relevant  expressions  for  the  perpendicular  wave  number  where  the  real  part 
of  the  frequency  is  equal  to  the  growth  rate  (i.e.,  kt)  and  the  perpendicular 
wave  number  appropriate  to  zero  growth  rate  (i.e.,  k$).  Hence,  kt  is  determined 
by  the  condition 

Re(wf)  =  Im(wf)  (Gl) 

while  k$  is  determined  by  the  condition 

0  *  Im(u»f)  .  .  (G2) 

Throughout  this  appendix,  the  dependence  of  on  k^  and  z  is  suppressed  for 
brevity  of  notation. 

In  the  strongly  driven  limit, 

(jj0,cu*  »  vj»vc  •  (G3) 

Equation  (G3)  is  most  likely  to  be  satisfied  in  early  or  intermediate  times 
following  a  high-altitude  nuclear  explosion.  The  combination  of  (3.36),  (3.37), 
and  (G3)  gives 


Re((A)r)  %  o 


Im(wf ) 


(•-?  4) 


(G4a) 


(G4b) 


«.* —  — V  --  -  * 


From  (Gl)  and  (G2),  it  follows  that 


.  ■  m 


(G5a) 


k$  -  (2)1/2kt 


(G5b) 


Hence,  ks  is  a  factor  of  (2)i/  larger  than  kt,  irrespective  of  plasma 
parameters. 

The  limit,  when  ion-neutral  collisional  effects  and  ion-viscous  effects 
are  significant,  is  important  for  late-time  phenomenology.  If 

vi  »  |ujf|  ,  |u)0|  ,  |ujt|  ,  vc  (G6) 

then  (3.36)  and  (3.37)  give 


Re(u4:) 


“o^i 

,  ,*  + 
ui  2, 

Vi 


(G7a) 


Im(w,) 


(G7b) 


The  first  and  second  terms  on  the  right-hand  side  of  (G7a)  require  the 
appropriate  consideration  of  ion-neutral  diffusion  and  ion  inertia, 
respectively.  The  first  term  on  the  right-hand  side  of  (G7b)  is  the  maximum 
growth  rate  in  the  collisional  limit  [cf.  (3.5)].  The  second  term  is  a  result  of 
the  stabilizing  role  played  by  ambipolar  diffusion. 

If  the  contribution  of  ion-viscosity  to  v..  is  small  relative  to  the 
contribution  of  ion-neutral  collisions,  i.e.. 


then 


kt  =  -Cj  *  ccf  ♦  t|)1/2 


( G9a ) 


ks  =  c2 


(G9b) 


In  (G9a) 


.  /  <dT.  fl.u. 
r  _  1  I  B  10  .  1  10. 

^  I  _  ^  I  P\  "r'  o  1 


1  2\mifliLiDc  Liv2n0c 
If  c2  «  then  (G9a)  reduces  to 


=  Vioy 

2  LivinDc 


( G10 ) 


kt  35  ks 


(Gil) 


Equation  (Gil)  implies  that  the  real  part  of  the  frequency  tends  to  be  small 

2  2 

compared  to  the  maximum  growth  rate,  If  then  (G9a) 

reduces  to 


k  =  _ _.c_ 

t  \>.  kdT . 
in  B  io 


(G12) 


This  value  for  kt  is  much  smaller  than  k$  and  indicates  that  the  real  part  of  the 
frequency  can  be  comparable  to  or  larger  than  the  maximum  growth  rate, 
tj0ou)i^ivin*  Ration  (G12)  is  precisely  what  would  be  obtained  in  a  static 
situation  if  an  ambipolar  electric  field  of  strength,  l9e|/vjn»  were  to  be 
balanced  by  an  ion-pressure  gradient.  Sperling  and  Glassman  [1983]  show  that 
(G12)  provides  a  better  estimate  for  kt  than  (Gil)  for  spread-F  and  barium 
clouds. 

If  the  major  contribution  to  the  effective  ion-collision  frequency  is 
from  ion  viscosity  rather  than  from  ion-neutral  collisions,  then 


vi  *  kiev  » 


(G13) 
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Sufficiently  large  values  for  kj^  are  required  for  (G13)  to  be  valid.  Frofr 
(G2)  and  (G7b)  it  then  follows  that 


/^iuiov 

|lM 


(614) 


If  vc  »  a)*,  then  (Gl)  and  (G7a-b)  give 


f  ft .  U  • 

i  loy 


,1/4 


(G15) 


The  expressions  for  k$  and  kt  in  (G14)  and  (GI5)  are  equal  because  the 
characteristic  growth  rate,  a>0ui*/^ ^ ,  is  much  larger  than  Re(w)  (scot).  If 
vc  «  co*,  then  (Gl)  and  (G7a-b)  give 


kt  = 


miai“loy 
IC8^io0v 


,1/3 


( G16 ) 


When  (G16)  is  appropriate,  ks  is  substantially  smaller  than  kt  and  Re(u>)  (soi) 
can  be  much  larger  than  the  characteristic  growth  rate,  <joQcot 

AiV 

In  general,  kt  %  k$  when  vc  »  cot  and  the  real  part  of  the  frequency  is 
unimportant.  However,  kt  «  ks  when  \>c  «  cot.  The  real  part  of  oo  then  gives  a 
significant  contribution  to  k^. 

The  various  approximate  expressions  described  in  this  appendix  (for  u>, 
k.,  and  k  )  are  summarized  in  Tables  Gl,  G2,  and  G3. 


asymptotic  expressions  for  u»f  as  deduced  from  (3.36 


Expression  I  Criteria  for  Validity 


Expression 
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NUMERICAL  INVESTIGATION  OF  THE  RELATIVE  IMPORTANCE 
OF  VARIOUS  ION  DRIFTS  ON  STRUCTURING 

To  evaluate  the  relative  importance  of  ion-Pedersen  drifts  and  the 
combination  of  ion-gravity  and  ion-curvature  drifts  on  possible  structuring 
during  late-times  following  a  high-altitude  nuclear  burst,  a  fifth  case  is 
specified  in  Table  3.2.  The  table  shows  that  Case  5  differs  from  Case  1  in  that 
Case  5  neglects  ion-gravity  and  ion-curvature  drifts  and  so  only  ion-Pedersen 
drifts  in  the  ionosphere  can  cause  structuring  instabilities.  Case  1  includes 
ion-gravity,  ion-curvature,  and  ion-Pedersen  drifts.  For  Case  1,  instability 
can  be  driven  at  high  as  well  as  low  altitude.  Of  the  two  cases.  Case  1  is 
clearly  more  realistic  than  Case  5. 

Figures  3.3a  and  Hla  are  plots  of  | | /k±  versus  2  (km)  for  Cases  1  and 
5,  respectively.  Curves  are  shown  in  the  figures  for  different  kj_  (km"'*').  As  in 
the  nonnuclear  examples  of  Sperling  and  Glassman  [1983],  the  striation  eigen- 
modes  become  more  localized  along  the  geomagnetic  field  as  the  perpendicular  wave 
number  becomes  larger.  Ion-polarization  currents  are  essential  to  the 
localization  at  higher  altitudes  where  ion-Pedersen  currents  are  negligible 
[Sperling,  1983e].  Evidently,  in  late-time  HANE  plasmas,  localized  disturbances 
along  the  geomagnetic  field  generally  lead  to  a  localized  response  and  so  the 
geomagnetic-field  lines  are  not  equi potentials.  The  major  qualitative 
difference  between  Figures  3.3a  and  Hla  is  the  tendency  for  the  maximum  value  of 
I  Ely | /k,  to  occur  at  lower  altitudes  in  Figure  Hla  and  at  higher  altitudes  in 
Figure  3.3a.  This  suggests  that  ion-Pedersen  drifts  are  not  always  the  primary 
drivers  of  instability  in  nuclear  plumes.  Rather,  the  ion-gravity  and  ion- 
curvature  drifts  at  altitudes  above  the  neutral  atmosphere  can  be  the  major 
drivers  of  instability.  This  is  verified  by  Figure  Hlb  which  is  a  plot  of 
a)  (s'1)  versus  (km"1)  for  both  Cases  1  and  5.  Figure  Hlb  shows  that  growth 
rates  for  Case  1  are  larger  than  those  for  Case  5  for  all  perpendicular  wave 
numbers  except  right  near  stability.  Indeed  the  growth  rates  for  Case  1  a*"? 
larger  than  the  maximum  possible  growth  rate  for  Case  5  (i.e.,  To  in  Table  3.2), 
for  1  km"1  <  k^  <  14  km"1. 
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Figure  Hla.  |Ely(/k  versus  z  (km)  for  various  k±  (km-1)  and 

Case  5.  Each  of  the  curves  are  in  normalized 
dimensionless  units. 


Figure  Hlb.  w(s-1)  versus  kx  (km-1)  for  Cases  1  and  5.  The  dashed 
line  shows  the  appropriate  eigenvalues  with  the  neglect 
of  electromagnetic  effects.  The  dotted  line  shows  the 
eigenvalues  for  Case  1  neglecting  viscosity. 


The  dashed  lines  in  Figure  Hlb  are  the  growth  rates  obtained  when 
electromagnetic  effects  are  neglected.  In  agreement  with  •  the  nonnuclear 
examples  of  Sperling  and  Glassman  [1983],  the  neglect  of  electromagnetic  effects 
results  in  reduced  growth  rates  for  smaller  values  of  kx.  Evidently,  electro¬ 
magnetic  effects  hinder  electrostatic  short-circuiting  of  the  instability 
[Sperling,  1983e].  When  electromagnetic  effects  are  important,  the  growth  rate 
is  a  relatively  weak  function  of  k^.  Indeed,  our  numerical  results  also  indicate 
that  lEjyl/^  is  a  weak  function  of  k^  when  electromagnetic  effects  are  important 
(see  Figures  3.3a  and  Hla). 

For  both  cases,  ion  viscosity  is  considered  and  contributes  to  the  curves 
indicated  by  the  solid  and  dashed  lines  in  Figure  Hlb.  For  Case  1,  eigenvalues 
have  also  been  calculated  neglecting  ion  viscosity.  The  results  are  indicated  by 
the  dotted  line  in  Figure  Hlb.  It  is  clear  that  ion  viscosity  plays  a  vital 
role  in  structuring  for  larger  values  of  k^.  Equations  (G14)  and  (G15)  in 
Appendix  G  are  the  relevant  asymptotic  expressions  for  k$  and  kt, 
respectively.  Figure  Hlb  verifies  that  kt  =  k$. 
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